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CRACKING  OF  A  HOMOGENEOUS  HALF  PLANE 
DUE  TO  SLIDING  CONTACT 

Serkan  Dag  and  Fazil  Erdogan 

ME-MECH.  Department,  Lehigh  University,  Bethlehem,  PA  18015 


Abstract 

In  this  report  the  initiation  and  subcritical  growth  of  surface  cracks  in  homogeneous 
materials  due  to  sliding  contact  are  considered.  The  load  is  applied  through  a  rigid  stamp  with  an 
arbitrary  profile.  The  problem  is  formulated  and  solved  under  the  assumptions  of  plane  elasticity 
and  Coulomb  friction.  The  stress  state  on  the  surface  of  the  half  plane  is  analyzed  in  the  absence 
of  any  cracks  in  order  to  determine  the  crack  initiation  force  and  angle.  A  surface  crack  is  then 
introduced  and  the  resulting  coupled  crack/contact  problem  is  investigated.  Extensive  results  are 
presented  regarding  mostly  the  mixed  mode  stress  intensity  factors  which  constitute  the  main 
driving  force  for  the  subcritically  growing  cracks. 

1.  Introduction 

In  this  report,  the  edge  crack  problem  in  a  homogenous  half-plane  due  to  sliding 
contact  will  be  considered.  The  problem  geometry  is  shown  in  Figure  1.  A  homogeneous 
elastic  half-plane  is  indented  by  a  rigid  stamp  of  arbitrary  profile.  The  contact  area 
extends  from  x  =  a  to  x  =  b  at  the  surface  and  the  half-plane  contains  an  edge  crack  of 
length  d.  The  crack  is  perpendicular  to  the  boundary  of  the  half-plane.  The  contact  is 
assumed  to  transfer  both  normal  and  tangential  forces.  A  normal  force  of  P  is  applied  by 
the  rigid  stamp  and  sliding  contact  is  assumed  at  the  interface.  Friction  coefficient  is 
denoted  by  rj,  hence  the  tangential  force  transferred  by  the  contact  is  equal  to  rjP.  p  is  the 
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shear  modulus  of  the  half-plane,  k  is  the  Kolosov's  constant  and  for  this  plane  strain 
problem,  it  is  equal  to  3  —  Au  where  v  is  the  Poisson's  ratio. 


In  this  report,  the  two  dimensional  elasticity  problem  will  be  formulated  and  solved  using 
Navier's  equations.  By  the  use  of  Fourier  Transforms,  the  problem  will  be  reduced  to  a 
system  of  three  integral  equations,  and  these  equations  will  then  be  solved  numerically  to 
determine  the  stress  intensity  factors  at  the  crack  tip  and  contact  stresses. 

2.  Formulation  of  the  problem 

First,  we  will  express  the  boundary  conditions  that  must  be  satisfied  in  this  problem. 
There  are  mixed  boundary  conditions  at  the  surface  y  =  0  and  at  the  crack  plane  x  —  0. 
At  the  surface  of  the  half-plane  y  =  0,  shear  and  normal  stresses  are  zero  outside  the 
contact  area,  in  the  contact  area  normal  displacement  component  v(x,  0)  is  known.  At  the 
crack  plane  x  =  0,  crack  faces  do  not  transfer  any  normal  or  shear  stress  and  they 

displace  relative  to  each  other  in  normal  and  tangential  directions,  but  outside  the  crack, 

2 


i.e.  if  y  <  d,  there  is  no  relative  tangential  or  normal  displacement  at  x  =  0  plane.  The 
conditions  at  the  crack  plane  are  valid  only  if  the  crack  is  open,  under  some  loading 
conditions  crack  faces  may  contact.  However,  the  problem  will  be  solved  by  assuming  an 
open  crack  and  a  negative  mode  I  stress  intensity  factor  will  indicate  that  crack  closure 
occurs  for  that  loading  condition.  There  are  also  additional  boundary  conditions, 
integration  of  the  normal  stresses  in  the  contact  area  must  give  the  total  force  P  applied  to 
the  contact  and  stresses  must  vanish  as  (x2  +  y2)  — >  oo.  Since,  sliding  contact  is 
assumed,  shear  stress  in  the  contact  area  is  given  by  Coulomb's  Law.  All  the  boundary 
conditions  can  be  expressed  mathematically  in  the  following  form, 


GyyiXi  0)  —  0, 

x  <  a  and  x  >  b, 

(la) 

axy(x,  0)  =  0, 

x  <  a  and  x  >  b, 

(lb) 

K+ldxv{x’0)=f{x) • 

a  <  x  <  b. 

(lc) 

0"xj/(x,  0)  Tj<7yy{Xi  0), 

a  <  x  <  b, 

(Id) 

Czx(0,  y)  =  0, 

d  <  y  <  0, 

(le) 

Oxy{ 0,  y )  =  0, 

d  <  y  <  0, 

(If) 

fb 

/  cryy(x,  0)dx  =  —  P, 

J  a 

(lg) 

(^yy?  & xx)  *  0  (*£ 

+  y2)  ->  oo, 

(lh) 

where  v(x,  y )  is  the  normal  displacement  component  and  /(x)  is  a  known  function.  Note 
that  in  (Id)  cryy{x,  0)  is  always  negative  and  if  the  friction  coefficient  is  taken  as  positive, 
the  shear  force  transferred  by  the  contact  will  be  directed  as  shown  in  Figure  1.  The 
problem  will  be  formulated  using  three  unknown  functions  and  using  superposition. 
Considering  Figure  1  and  boundary  conditions  (1),  following  unknown  functions  can  be 
defined: 
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fV  (u(0+,y)  u(0  ,y))-fi{y), 

K  +  loy 

d  <  y  <  0, 

(2a) 

q  ( v(0+,y )  v(0  ,y))  =  /2(y), 

k  + 1  oy 

d  <  y  <  0, 

(2b) 

(«£?  0) 

a  <  x  <  b, 

(2c) 

where  u  and  v  are  the  displacement  components  in  x  and  y  directions  respectively.  The 
scheme  that  will  be  employed  in  the  formulation  is  shown  in  Figure  2.  Problem  1  is  the 
contact  problem  without  crack  and  in  this  case  the  stresses  and  displacements  will  be 
obtained  in  terms  of  the  unknown  contact  stress  which  is  given  by  (2c).  In  problem  2 
stress  and  displacement  fields  will  be  obtained  in  terms  of  the  relative  displacement 
derivatives  of  the  crack  faces  which  are  given  by  (2a)  and  (2b).  This  can  also  be  viewed 
as  the  distributed  dislocation  solution  for  the  homogenous  half-plane.  Note  that  f\  and  / 2 
are  zero  in  problem  1  and  fz  is  zero  in  problem  2.  The  total  stress  and  displacement  fields 
for  the  original  problem  can  then  be  obtained  by  summing  the  solutions  of  the  problems  1 
and  2,  and  satisfying  the  boundary  conditions  of  the  original  coupled  problem.  The 
stresses  and  displacements  may  then  be  expressed  as 


y)  =  (x,  y)  +  e-fix,  y). 

i,j  =  x  or  y, 

(3a) 

u(x,  y)  =  u(1)(x,  y)  +  u(2)(x,  y), 

(3b) 

v(x,  y)  =  u(1)(x,  y)  +  x(2)(x,  y). 

(3c) 

Superscripts  (1)  and  (2)  stand  for  problems  (1)  and  (2)  respectively.  The  stress  and 
displacement  fields  given  by  (3)  must  satisfy  the  conditions  given  by  (1).  In  the  following 
sections  problems  1  and  2  will  be  formulated  to  derive  the  stress  and  displacement  fields 
in  terms  of  a  certain  set  of  unknown  functions. 
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Figure  2:  Superposition 


2.1  Problem  1  -  The  Contact  Problem 


Geometry  of  the  contact  problem  is  shown  in  Figure  3.  In  this  section  stresses  and 
displacements  will  be  obtained  in  terms  of  the  unknown  contact  stress,  crj,y(a;,0)  at  the 
surface.  The  equations  of  equilibrium  can  be  written  in  the  following  form, 


d°xx  d°xy  _  0 
dx  dy 


dy  dx 


=  0. 


(4a) 

(4b) 


Assuming  plane  stress  or  plane  strain  and  small  deformations,  for  the  linear  elastic 
medium  considered  Hooke's  law  becomes 
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a„(x,!,)  =  ^T{(K  +  l)^  +  (3-K)|j},  (5b) 

+  <*> 

Contact  problem  is  a  plane  strain  problem,  hence  for  this  case  k  is  equal  to  3  -  iu. 
Substituting  (5)  in  (4a)  and  (4b),  governing  equations  for  the  displacements  can  be 
obtained  as  follows, 


d^u  <92u  <92u 

(*+1fe  +  <*-V  +  2^  =  0' 


.  -,\®2v  t  o  d2u  n 

^K~^d^  +  ^  +  1^df+2dxdy  ~  ' 


(6a) 

(6b) 


Considering  Fourier  Transformation  in  x,  u  and  v  can  be  expressed  in  the  following  form, 


1  f00 

u(x.y)  =  —  U  (w,  y)exp(iux)du, 

2tt  J — oo 


(7a) 


6 


(7b) 


1  [°° 

v{x,  y)  =  —  /  v(w>  y)exp(iwx)da;, 

2t rj-co 

where,  i  =  \/  —  l  and  U  (lo,  y)  and  V (a»,y)  are  Fourier  transforms  of  u(x,  y)  and  v(x,  y) 
in  x,  respectively.  Substituting  (7)  in  (6)  following  ordinary  differential  equations  are 
obtained 


d2TJ  dV 

-  u2(k  +  1)U  +  (k-  1)-^" +  2iu~d^  =  °’ 

—  u?(k  -  1)V  +  (k  +  l)-T—r-  +  2 iu—  =  0. 

dy 2  dy 


(8a) 

(8b) 


Solving  equations  (8),  substituting  the  solution  in  (7)  and  considering  the  regularity 
condition  given  by  (lh),  displacements  can  be  expressed  as  follows: 


1  f00 

u(x,y)  =  tt  {MC\  +  A2C2y)e^p{\to\y  +  iux)du, 

2  nj-oo 

1  f°° 

v{x,  y)  =  —  /  (Cl  +  C2y)exp(|u;|y  +  itoxjdu. 

2tc  J  _00 

Note  that  in  (9)  y  <  0,C\  and  C2  are  unknown  constants  and  A\  and  A2  are  given  by 


(9a) 

(9b) 


^  zk  C2  +  i\u>\ 


oj  C\  u 


Ao  = 


JW 


u 


(10a, b) 


Substituting  (9)  in  (5)  stresses  can  also  be  obtained 

axx(x, y)  =  -  7^  J  [2|w|Ci  +  ((«  +  3)  +  2|w|y)C2]  exp(|u;|y  +  iux)du, 

ayy(x,y)  =  7^J  [2|w|Ci  +  ((«  -  1)  +  2|u;|y)C2  exp(|cu|y  +  iux)du, 

axy{x,y)  =  ~  [  \2uCi  +  ((«  +  1)M  +  2u2y)—]exp(\u\y  +  iux)du. 
27T  J-m  l  u  J 


(11a) 

(lib) 

(11c) 
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The  constants  C\  and  C2  will  be  determined  using  the  conditions  (la),  (lb)  and  (Id)  and 
stresses  and  displacements  will  be  expressed  in  terms  of  the  contact  stress  ayy(x,  0). 
Considering  boundary  conditions  (la),  (lb),  (Id)  and  (2c)  we  can  write 


ayy{x,  0 )  =  7^  J  2MCi  +  (*  “  1)^2]  exp (iux)du  =  |  ^ 


f3(x),  a  <  x  <  b 
0,  x  <  a,  x  >  b 


axy(x,  0)  =  ^  f  [2 uCi  +  (k  +  1 ) <-^2]  exp ( iux)du> 


rjh(x),  a  <  x  <  b 
0,  x  <  a,  x  >  b 


(12b) 


Taking  Fourier  transforms  of  both  sides  in  (12)  we  obtain  the  following  equations, 

fj,(2\u>\Ci  +  (k-  1  )C2)  =  f  fz(t)exp(  -  iut)dt, 

J  a 

jx^2iojCi  +  »(k  4-  1)MC2)  =  J  77/3  (t)  exp  (  -  iut)dt. 


(13b) 


Ci  and  C2  are  determined  in  terms  of  /3  solving  the  linear  system  of  equations  (13)  and 
they  are  substituted  in  (11)  to  derive  the  expressions  for  the  stresses.  These  operations  are 
carried  out  using  the  symbolic  manipulator  MAPLE,  some  integrals  are  evaluated  in 
closed  form  and  following  results  are  obtained: 

,  x  2  fbr]{t-x)z -y{t-xf  nA^ 

y)  =  "J.  ’  (  > 


<7yy(X 


>y)  =  l[ 

^  J  a 


/(*, y)  =  lf 

7T  Ja 


«Ja  (y2  +  (t  —  x)2) 

2  fb  r]y2(t  -x)-y3 

2  fb  -  rjy(t  -  xf  +  y2(t  -  x) 


(14b) 


(y2  +  (t-x)2Y 


Mt)dt. 
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the  resulting  system  of  ordinary  differential  equations  the  expressions  for  stresses  and 
displacements  for  both  half-planes  x  >  0  and  x  <  0  can  be  obtained  by  following  the 
procedure  described  below. 


Half-Plane,  x  >  0 


In  the  following  equations  superscript  4  stands  for  the  half-plane  x  >  0. 


-oo  L 
poo 


t*(4)(*,  y)  =  ir  —°l4)  +(-  +  x )  a24)  exP(  -  M*  +  ™y)d“>  (15a) 

2?r  J-oo  a>  \u  u  J 


v{4){x,  y)  =  +  C2(4)x)ex p(  -  \u\x  +  iojy)du, 

/oo  -  (7^4)  -j 

2 ujC[A)  +  ({k  +  1)M  +  2<Jx)^L-  exp(  -  \u\x  +  iuy)du, 

■oo  *-  ^  ^ 


(15b) 


•  /*00  r  /^(4)  -| 

$  (X,  y)  =  7^J  [2uCi4)  +  ( -  (3  -  k)  M  +  2w2a:)  exP( 


(15c) 

cu|x  +  iuy)du, 

(15d) 

/OO  “i 

[2|c^|C'1(4)  +  ((«  -  1)  +  2|u;|x)C'2(4)  exp(  -  |u/|®  +  twy)dw, 

C15c) 

where  and  C^A>  are  unknown  constants.  Similarly  the  solution  for  x  <  0  is  obtained 
in  the  following  form: 


Half-Plane,  x  <  0 


In  the  following  equations  superscript  3  stands  for  the  half-plane  x  <  0. 


v^,y)  = 


U  \  U  U  ' 


exp(\u>\x  +  iuy)du,  (16a) 
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Figure  5:  Superposition  for  the  crack  problem 

where  C®  and  are  unknown  constants.  The  boundary  conditions  for  the  infinite 
plane  problem  are  in  the  following  form. 
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o**(0  ,y)  =  4x(0,r/), 

—  oo  <  y  <  oo, 

(17a) 

<4?(o,y)  =  c4?(o,y), 

—  oo  <  y  <  oo, 

(17b) 

u(3)(0,y))  =  < 

f  /i(y)» 
10, 

d  <  y  <  0 
y  <d,y>  0 

(17c) 

u(3)(0,?/))  =  <j 

7a(y). 

d  <  y  <  0 
y  <d,y>  0 

(17d) 

Using  the  expressions  (15),  (16)  and  boundary  conditions  (17),  the  four  unknown 
constants  are  determined  in  terms  of  fi  and  f2  using  the  symbolic  manipulator  MAPLE. 
The  constants  are  then  substituted  in  the  expressions  for  stresses  and  displacements,  and 
the  stresses  and  normal  displacement  derivative  for  the  infinite  plane  containing  a  crack 
are  obtained  as  follows: 


°>’y)=U 


i  [°(y-t){x2  -  (y-t)2) 


71  Jd  (x2  +  (y-t)2)2 

1  f°x(x2  +  3 (y  -  tf ) 


+  i/  7  o  V 

KJd  (x2  +  (y-t)  ) 


<Jxi(x,y)  = - 

7 r 


1  f°{y-t){3x2  +  (y-t)2) 


* Jd  (x2  +  (y  —  t)2)2 
1  f°x(-x2  +  {y-t)2) 


irJd  (x2  +  (y-t)2)2 
1  f°x(x2  -{y-  t )2) 


f2(t)dt, 


^,y)  =  i  ;a  . 

nJd  (x2  +  (t/-f)) 


+1  f0(y-^)(^2  -  (s/-*)2) 

nJd  (x2  +  (y  -  t)2)2 


f2{t)dt, 


(18a) 


(18b) 


(18c) 
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(1 8d) 


ivii){x’y)  =  ^-Jd 


—  f 

4ar  pJd 


'  (1  -  k)x 3  -  (3  +  K)x(y  -  tf 


+ 


( x 2  +  (y- 1 )2)2 
{y  -  t )  [(1  -/t):r2-(3  +  K)(y-t)2] 


2\2 


(*2  +  (2/  -  t)2) 


where  superscript  (i)  stands  for  the  infinite  plane.  Half-plane  y  <  0  solution,  can  now  be 
superimposed  on  the  solution  for  the  crack  in  the  infinite  plane  to  obtain  the  required 
solution.  The  half  plane  y  <  0  solution  is  given  by  equations  (9)  and  (1 1)  in  terms  of  two 
unknown  constants  C\  and  C2.  These  constants  are  determined  in  this  case  by  imposing 
the  following  free  surface  conditions: 


(Jyy(X,  0)  —  0, 

—  OO  <  X  <  oo, 

(19a) 

aXy(x,  0)  =  0, 

—  OO  <  X  <  oo. 

(19b) 

Note  that  these  conditions  must  be  satisfied,  by  the  total  solution,  i.e.,  the  sum  of  the 
solutions  for  infinite  and  half-planes  have  to  satisfy  the  boundary  conditions.  The 
constants  C\  and  C2  are  determined  using  the  symbolic  manipulator  MAPLE,  by 
imposing  the  above  mentioned  boundary  conditions,  the  stresses  and  normal 
displacement  derivative  for  a  crack  in  a  half  plane  as  depicted  in  Figure  4,  are  then 
obtained  as  follows: 

oxx(x,y)  =  ^  Sxxl(x,y,t)fi{t)dt  +  Sxx2(x,y,t)f2(t)dt  j  +aQ(x,y), 

(20a) 

ayy(x,y)  =  Syyi(x,y,t)fi(t)dt  +  Syy2{x,y,t)f2(t)dt j  +a^(x,y), 

(20b) 

axy{x,y)  =  i  j  Sxyi(x,y,t)f1(t)dt  +  Sxy2(x,y,t)f2{t)dt^  +a^(x,y), 

(20c) 
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Mm)  .  _L {  j\l(x, y,t)hm  +  [sax, v, t)h(t)dt }  + 

(20d) 


Sxxi(x,  y,  t )  -  {l6yi(y  +  tf  -  2(y  +  t)(y2  +  10 yt  +  3 t2)  ((y  +  tf  +  x2) 

+  (3y  +  5t)((y  + 1)2  +  x2)2|/((y  +  t)2  +  x2)  ,  (20e) 

SXx2{x,y,t)  =  [l6xyt(y  +  tf  -  2x(3t2  +  Ayt  -  y2)((y  +  tf  +  x2) 

-x((y  +  tf  +  x2f}/((y  +  tf  +  x2)3,  (2  Of) 

Syy i(x,  y,i)  =  {  -  16yt(y  + 1 f  -  2(y  + 1)(  -  y2  -  6yt  + 12)  ((y  +  tf  +  x2) 

+  (  -  y  + 1)  ((y  +  tf  +  x2)2}/((y  + 1)2  +  x2)3,  (20g) 

Syy2(x,  y,  t)  =  {  -  16xy*(y  +  tf  -  2 x(y2  +  t2)((y  +  tf+  x2) 

-  x((y  +  tf  +  x2)2}/((y  +  tf  +  x2)3,  (20h) 

SXyi(x, y, t)  =  {  -  16xyt(y  +  tf  +  2x(4yt  +  y2  + 1 2)  ((y  +  tf  +  x2) 

-x((y  +  tf+x2f}/((y  +  tf+x2f,  (20i) 

Sxy2(x,  y,  t)  =  {l6y*(y  +  tf  -2 (y  +  t)(-y2+  6 yt  + 1 2)  ((y  + 1)2  +  x2) 

+  (*  -  y)((y  + *)2  +  z2)2}/((y  +  *)2  +  *2)3>  (2°j) 

^i(x,y,  t)  =  x|(k  +  3)(y4  —  8yt3)  —  6(k  +  8)y2t2  +  (k  —  l)x4 
+  2(k  +  l)x2(y2  -  t2)  -  (1  +  3 K,)t4 

+  8 yt(x2  -  2y2)|/((y  +  tf  +  x2)3,  (20k) 
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•^(z,  *)  =  {(*  +  %5  +  2(K  +  ^(y3  +  *3)  +  («  -  l)^4(y  -  0 

+  6(«  -  3 )yxH(y  +  t)  +  2(9 k  +  23 )t2yz  +  2(11«  +  21  )t3y2 
+  (13k  +  5  )tAy  +  (7k  +  21)fy4  +  (1  +  3  K)t®}/((y  +  tf  +  x2)3.  (201) 


2.3  Solution  for  the  coupled  crack  -  contact  problem 


Since  the  expressions  for  the  contact  and  crack  problems  are  derived,  stresses  and 
normal  displacement  derivative  for  the  coupled  crack  and  contact  problem  can  now  be 
obtained  by  using  equations  (14)  and  (20).  For  the  coupled  problem  stresses  on  the  crack 
plane,  and  normal  displacement  derivative  at  the  contact  surface  are  found  to  be 
'  i  i  2 1  At2 


o. 


+ 


+  y  ( t  +  y )  (t  +  y )' 


+ 


2  fb  yt3 
Kja  ( y 2 


w  -  yt 


(y2  +  t2)2 

0/1 


oo  <  y  <  0, 


+ 


2 1 


4 12 


y  t  +  y  (t  +  y)  (t  +  y)' 


■)/2(t)dt 


+ 


2  rb 

KJa 


-  rjytf  +  ylt 

(: y 2  + 1 2)2 


fz(t)dt. 


4f*  8  ,  4  f°  (2x  4  f 


—  oo  <  j/  <  0, 
3  t3 


(x2+t2)2yiW  ‘  7T Jd  (x2+t2)2 
k-1^  N  ,  1  /-6/3(i)  ^ 

-  ’'rrr/3(x)  +  *L  —xdt’ 


f2(t)dt 


OO  <  X  <  oo. 


(21a) 


(21b) 


(21c) 


15 


3.  Integral  equations  and  singular  behavior  of  the  unknown  functions 


Using  the  boundary  conditions  (lc),  (le)  and  (If)  and  equation  (21),  following  system 
of  integral  equations  can  be  obtained  to  determine  the  unknown  functions  f\ ,  f2  and  fz 
for  the  coupled  crack  and  contact  problem, 

<7M(0 ,y)  =  -  f^dt  +  -  [°Ku(t,y)fi(t)dt  +  -  f  Kxz(t,y)h(t)dt  =  0, 

7T  Jd  t-y  TT  Jd  IT  Ja 

d  <  y  <  0,  (22a) 

Cxj^O,  y)  =  -  f  ~f2 -- dt  +  -  (  K2i{t,y)f2{t)dt  +  -  f  K2z{t,y)fz(t)dt  —  0, 

TT  J  d  t  y  IT  J  £  TT  J  a 

d<y<  0,  (22b) 

!-»(*, 0)  =  -  f  K3l(t,x)h(t)dt  +-  f  Kn{t,x)h(t)dt  -  n^-h(x) 

(«  +  1)  dx  rcjd  7 Tjd  K+l 

+  -  t^-dt  =  f(x),  a  <  x  <  b.  (22c) 

KJa  t  ~  X 

The  kernels  are  given  by 


r,  ,  ,  r,  ,  ,  1  2 1 

Ku(t,y)  =  K22(t,y)  —  — - H 


4r 


t  +  y  ( t  +  y Y  (t  +  y) 


3  ’ 


Kis(t,y)  = 
K2z(t,y)  = 


2(r)t3  -  yt2) 

(t/2+t2)2  5 

2(  -  rjyt2  +  y2t) 


Kzi(t,x)  =  - 
K32(t,x) 


{y2  + 1 2) 

4  t2x 

(x2+t2)2’ 
4 13 


2  ’ 


(x2+f2)2’ 


(23a) 

(23b) 

(23c) 

(23d) 

(23  e) 
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All  the  singular  integral  equations  contain  a  Cauchy  integral  and  also  the  kernels  Kn, 
K22i,  A13  and  K23  become  singular  as  t  and  y  simultaneously  go  to  zero.  Similarly,  the 
kernels  KZ\  and  K32  become  singular  as  x  and  t  simultaneously  go  to  zero.  Since, 
axx{0 ,y),  0xy(O,y)  and  dv(x,0)/dx  are  bounded  at  all  points  in  their  respective 
intervals,  the  unknown  functions  may  have  a  special  behavior  in  the  form  of  a  power 
singularity  at  the  end  points.  In  this  section  the  singular  behavior  of  the  unknown 
functions  will  be  examined  using  a  function-theoretic  method.  In  this  analysis,  there  are 
two  cases.  If  b  <  0,  in  addition  to  the  Cauchy  singularities  only  terms  that  can  become 
singular  are  Ku  and  K22.  Other  kernels  are  bounded  at  all  points  of  their  respective 
intervals.  If  b  =  0,  all  of  the  kernels  have  to  be  examined  to  determine  the  singular 
behavior  of  unknown  functions.  We  will  begin  with  the  case  of  b  <  0  in  the  next  section. 

3.1.  Singular  behavior  of  the  unknown  functions  for,  b  <  0 

We  will  examine  the  possibility  of  a  power  singularity  in  the  unknown  functions  for 
the  case  b  <  0.  Initially  we  assume  the  following  form  for  the  unknown  functions, 

h(t)  =  Fi(f)(-f)ai(t-d)71,  d<t<  0,  (24a) 

hit)  =  F2(t)(  -  -  df\  d<t<  0,  (24b) 

h{t)  =  F3{t)(b  -  t)w(t  -  af,  a  <  t  <  b.  (24c) 

where  Fft),  (j  =  1, 2, 3)  is  Holder-continuous  in  its  respective  interval  and  it  is  assumed 
that  -  1  <  3fJ(o;i,a2, 71,72,  v,/3)  <  0.  Consider  now  the  following  sectionally 
holomorphic  functions 

=  1  (25a) 

TtJd  t~Z 
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Mz)  = ; 

7T  J  d  t  z 

J a 


t  —  Z 


dt . 


(25b) 

(25c) 


The  singular  behavior  of  ipj(z),  ( j  =  1, 2, 3)  around  the  end  points  is  in  the  following 
form,  (see,  for  example  Erdogan  [1]  or  Muskhelishvili  [2]), 


M*)=  -Fi(d)(-d) 


<*1 


exp(  —  7ri7i) 
sin(7T7i) 


1 


(z-^  +  Fl(0)(-dr-^+F-W, 


(26a) 


exp(  -  7T272) 


m*)  =  -  wx  -  dr  - rf)72 + ^(o)(  -  d) 


\12 


ra2 


sin(7T72) 


sin(7TO!2) 


+  F*(z), 


(26b) 


=  ~  Fz{a)(b  -  a)' 


;exp(-7ri/5)/_  p  ,  ^^(2-6)c 


sin(7r/?) 


-(z-af  +  F3(6)(6-ay 


sin(mu) 


+  *»• 


(26c) 


The  function  F*(z),(n  =  1,2,3)  is  bounded  everywhere  except  possibly  at  the  end 
points,  where  it  may  have  a  weaker  singularity.  Using  Plemelj  formulas  and  equations 
(26)  the  singular  behavior  of  the  Cauchy  integrals  at  the  end  points,  for  the  given  form  of 
unknown  functions  is  obtained  as  follows: 

I  f°iMdt  =  _  Jp1(d)(  -  d)Qlcot(7r7l)(Z/  -  d )71  +  Fi(0)(  -  d)7lcot(7ra1)(  -  y)ai 
nJd  t-y  “ 

+  Hi(y),  d<y<  0,  (27a) 

I  rM>  di  =  _  p  (d)(  -  d)Q2cot(7r72)(2/  -  d)72  +  F2(0)(  -  d)72cot(7ra2)(  -  2/)°2 
nJd  t-y 

+  H2{y),  d  <  y  <  0,  (27b) 


-  f  ^3-—df  =  -  Fs(a)(fe  -  a)wcot(7r/5)(x  -  a)^  4-  Fz(b)(b  -  a)^cot(7ru;)(&  -  x)w 

7T Ja  t-X 


+  Hz{x),  a  <  x  <b. 


(27c) 
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Hi(y),  H2(y )  and  H3(x)  are  bounded  in  their  respective  intervals,  and  at  the  end  points 
their  behavior  is  similar  to  that  of  F*(z )  in  equations  (26).  Now,  consider  equations  (25) 
and  assume  that  complex  variables  z\,  z2  and  23  satisfy  the  following  conditions: 


*i  i  (d  <  y  <  0), 
z2£{d<y  <  0), 
z3  £  (a  <  x  <  b). 


(28a) 

(28b) 

(28c) 


Under  these  conditions,  ^1(21),  ^2 (22)  and  ipz{zz)  are  holomorphic. 

1  f°  fi(t) 


VdOi)  =  -  f 

J  d 


t  —  Z\ 

1  r°  m) 


dt, 


Mz2)  =  l[  J~~dt, 

nJd  t-  z2 

Mzs)  =  lf  Mt)  •" 

^  J  a 


t-Zz 


Thus,  we  can  write, 
(29a) 

(29b) 

(29c) 


For  the  case  b  <  0,  other  than  the  Cauchy  integrals  whose  singular  behavior  are  given  by 
(27),  the  kernels  that  can  become  singular  are  K\\  and  K22  and  all  the  other  kernels  of  the 
integral  equation  remain  bounded.  It  can  be  seen  that,  as  y  and  t  simultaneously  go  to 
zero  Kn  and  K22  become  singular.  Equations  (29)  will  be  used  to  determine  the  singular 
behavior  of  Ku  and  K22  in  the  following  analysis.  Kn  can  also  be  expressed  in  the 
following  form: 


Kn(t,y) 


1  6  y _ 4y2 

(t  +  y)  +  (t  +  y)2  ( t  +  y )3‘ 


(30) 


Also  note  that, 


1 

(t  +  V? 


d 

dy 


(31a) 
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Then,  using  (29a)  and  (30)-(31)  we  can  write 


U  (jfi 

\fdKn(t,y)h{t)dt  =  -  M  -  y)  -  ~v)~  ~  v)-  (32) 


Using  (26a),  singular  behavior  of  this  term  near  y  =  0  can  be  expressed  as  follows: 

lfdKn(t,ymm  « -  -  <r<  -  »)**«». 


Similarly  for  K2 2  we  obtain, 

l  ^  -  »>"*(0)-  <33b) 


Note  that  instead  of  showing  the  bounded  terms,  s*  sign  is  used  in  (33).  Since  singular 
behavior  of  all  the  terms  are  determined,  using  (27)  and  (33)  the  singular  behavior  of  the 
integral  equations  given  by  (22)  can  be  expressed  as  follows: 


axx(0,  y)  a  Fi(0)(  -  d )71  |cot(m*i)  - 


2  af  +  4ai  +  1 

sin(7ro:i) 


(~y) 


Oil 


-  Fi(d)(  -  d)aicot(7T7i)(y  -  d)7\ 


(34a) 


(Jxy{0,y)  S£  1*2(0)  (  -  d) 


72 


cot(7ra2)  — 


2q2  +  4Q2  +  1 

sin(7ro;2) 


]  ( -  v)“ 


-  F2(d)(  -  d)a2cot(7T72)(y  -  d)72 , 


(34b) 


-^——v(x,  0)  s  -  Fz(a)(b  -  a)ucot(n/3)(x  -  af  +  F3(b)(b  -  a/cot(mu)(6  -  x)u 
K  +  \dx 

-r?~X7 Fz(x)(b  -  x)u (x  -  a)p .  (34c) 

K  +  1 


Now,  the  characteristic  equations  can  be  derived  by  using  (34)  to  determine  the  unknown 

exponents.  Multiplying  (34a)  by  (  —  y)~ai  and  letting  y  — *  0  we  obtain, 
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(35) 


2  aj  +  4a i  +  1 

cot(7rai) - .  - t —  =  0, 

sm(7rai) 

and  multiplying  (34a)  by  ( y  —  d)~11  and  letting  y  —*  d  we  have, 


COt(7T7i)  =  0. 


(36) 


It  is  assumed  that  Fi(0)  and  Fx(d)  are  not  equal  to  zero.  Similarly  using  equation  (34b) 
we  obtain 


cot(7ra2) 


2al  +  4a2  +  1 
sin(7ra2) 


=  0, 


(37a) 


COt(7T72)  =  0. 


(37b) 


Again  F2(0)  and  F2(d)  are  assumed  not  to  be  zero.  Multiplying  (34c)  by  (x  -  a)  0  and 
letting  i->awe  have  the  following  equation 

cot(7 r/3)  =  —  rj - (38) 

V  K  +  l 


and  multiplying  (34c)  by  (b  -  x)  w  and  letting  x  b  we  have  another  characteristic 
equation  to  determine  u 


cot  (mu)  =  r] 


k  —  1 

K+l 


(39) 


Fz(a)  and  Fs(b)  are  assumed  not  to  be  zero.  Equations  (35)  and  (37a)  have  no  roots  in 
the  interval  —  1  <  a\,  a2  <  0,  hence  the  unknowns  /i  and  /2  have  no  power  singularity 
at  the  end  point  at  y  =  0.  This  is  the  well  known  result  for  an  edge  crack  and  from 
equations  (36)  and  (37b)  we  obtain,  71  =  72  =  —  1/2.  Equations  (38)  and  (39)  give  the 
strength  of  singularity  at  the  ends  of  the  contact  area.  These  exponents  are  less  than  zero 
if  the  half-plane  is  in  contact  with  a  flat  stamp  as  shown  in  Figure  1.  If  the  half-plane  is 
indented  by  a  circular  stamp  at  the  end  points  contact  stress  is  zero  and  the  exponents  are 
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greater  then  zero.  Since,  there  is  no  singularity  at  y  —  0  for  f\  and  f2,  they  can  be 
expressed  in  the  following  form 


m 


0 t-d )1/2’ 


f2(t)  = 


m) 

(t-df2' 


(40a, b) 


In  this  case  the  asymptotic  behavior  of  the  sectionally  holomorphic  functions  ipi(z)  and 
il>2(z)  near  the  end  point  0,  can  be  written  as  follows  (see,  for  example  Muskhelishvili, 
[2]), 

■01  (z)  s  -P=ln(z),  ^2(2)  =  — y==ln(2)-  (41a, b) 

7 ry  —  d  Tty  —  a 


Singular  behavior  of  the  Cauchy  integrals  at  the  end  point  0,  then  can  be  written  as 


t-y 


t-y 


£U Q) 

7T  y/—d 

Fm 

7 ry  —  d 


M  -  y)> 
ln(  -  y). 


(42a) 

(42b) 


For  this  case,  integral  equations  must  be  examined  to  search  for  a  logarithmic  singularity. 
Considering  equation  (30),  the  only  term  that  causes  a  logarithmic  singularity  is  the  first 
term  and  we  can  write, 

_  I  f!Mdt  a  _  J>£Lin(  -  »).  (43) 

^  Jd  t  +  y  7T\/  —  d 


Hence,  in  integral  equation  (22a)  the  only  terms  that  have  logarithmic  singularity  are 
given  by  (42a)  and  (43),  but  as  can  be  seen  these  terms  cancel  each  other.  So  there  is  no 
additional  condition  to  eliminate  the  existence  of  a  logarithmic  singularity  in  (22a)  and 
similarly  it  can  be  shown  that  this  is  also  the  case  for  (22b). 
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3.2.  Singular  behavior  of  the  unknown  functions  for,  6  =  0. 


In  this  section,  first  we  will  examine  the  possibility  of  a  power  singularity  in  the 
unknown  functions  for  b  =  0  using  the  function-theoretic  method.  If  b  =  0,  ifn,  if 22, 
if13  and  if 23  become  singular  as  y  and  t  go  to  0  and  if 31  and  if 32  become  singular  as  x 
and  t  approach  0.  Hence,  all  the  kernels  must  be  examined  to  determine  the  strength  of 
singularity  at  the  point,  x  =  0,  y  =  0.  The  unknown  functions  can  expressed  in  the 
following  form, 


m  =  Gft)(-tnt-df\  d< 

f2(t)  =  G2(t)(-tnt-dy\  d< 

fs(t)  =  Gz (t)(  —  t)a{t  —  a)P ,  a  < 


t  <  0, 

(44a) 

t  <  0, 

(44b) 

t  <  0. 

(44c) 

The  function  Gj(t),  ( j  =  1, 2, 3)  is  Holder-continuous  in  its  respective  interval,  and  it  is 
assumed  that  -  1  <  §R(q:,7i  ,72,/?)  <  O.The  definition  of  the  sectionally  holomorphic 
functions  is  given  by  equations  (25)  and  the  singular  behavior  of  these  functions  near  the 
end  points  is  given  by, 


~  + Gi<0)(  -  + G:w’ 


(45a) 


=  -  G’M<  -  ^  +  G2(0)(  -  O'sSsj  ■ +  G^' 


Ip3{z)  =  -  G3(o)(  -  a) 


•%  -  of  +  G»(0)(  -  of  -|-£-r  +  GHz). 

Sin(7 xp)  Sin(7TQ;) 


(45b) 


(45c) 
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G'n(z)  is  bounded  everywhere  except  possibly  at  the  end  points  where  it  may  have  a 
weaker  singularity.  Using  Plemelj  formulas  and  equations  (45)  singular  behavior  of  the 
Cauchy  integrals  at  the  end  points  can  be  expressed  as, 


-  f  dt  =  -  Gj(d)(  -  d)acot(7T7i)(y  -  d )71  +  Gi(0)(  -  d)7lcot(7ro;)(  - 

njd  t-y 

+  d  <  y  <  0, 

-  [°ML idt  =  -  G2(d)(  -  d)acot{nl2)(y  -  d)72  +  G2(0)(  -  d)72cot(7ra)(  - 
7 rjd  t-y 

+  H2(y),  d  <  y  <  0, 


vT 


(46a) 


vT 


(46b) 


-  f  hQl dt  =  -  G3(a)(  -  a)acot(7r/3)(x  -  af  4-  Gs(0)(  -  a)^cot(7ra)(&  -  x)a 

TTJa  t-X 

+  H3(x),  a  <  x  <  0.  (46c) 


If  we  assume  complex  variables  z\ ,  z2  and  z2  satisfy  the  following  conditions 


z\  i  (d  <  y  <  0), 

z2£(d  <y  <  0), 
z2  (a  <  x  <  0), 


(47a) 

(47b) 

(47c) 


then,  ipi(zi),  ip2(z2)  and  ^3(23)  are  holomorphic  and  equations  (29)  are  valid.  In  order  to 
determine  the  asymptotic  behavior  of  the  integral  equations,  the  kernels  will  be  expanded 
into  partial  fractions  in  the  following  form, 


Kn(t,y)  =  K22(t,  y)  = 

i 

2  (t  +  iy) 


_  1  6  y _ 4  y2 

(t  +  y)  (t  +  y)2  (t  +  yf  ’ 

1  jy _ iy  \ 

t  -  iy  +  2(t  -  iyf  2(t  +  iyf' 

j _ y _ y 

2(t  —  iy)  2(t  +  iy)2  2  (t-iy)2’ 


(48a) 


(48b) 
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#23  (t,y)  —  i)[  —  ~  —  + 


y 


y 


+ 


2(f +  iy)  '  2 (t-iy)  2{t  +  iyf  2 (t-iyf 

iy  iy 


2 (t  +  iy)2  2(i  -  iy) 


#3i(t,x)  = 

#32  (*,  ®)  = 


—  2 


X  2 

+ 


X 


(f  +  ari)  (i  +  xi)2  '  (i-xi)2’ 

2 


ix  2  ix 

-  +  -  :  + 


t  +  xi  (t  +  xi)  t-xi  (t  -  xi) 


.•\2  ‘ 


(48c) 

(48d) 

(48e) 


(49a) 


Using  (29)  and  (48)  we  can  write  the  following  expressions 

^  J  Ku(t,y)fi(t)dt  =  --0i(-y) -6y^^!(-y)-2y2^2/>i(-y). 

^  J  K22(t, y)f2(t)dt  =  -i>2 ( - y)  -  -v)~  2y2 ~  y ')»  (49b) 

^  J  Ki3(t,  y)fz(t)dt  =  r\ (rp3(  -  iy)  +  ipz{iy)  +  +  VhC  -  iy 


(49c) 


+  ^Mw)  -  M  ~  *y))  ~  ~  ^  ~ iy '))» 

^ J  K2z{t,y)h{t)dt  =  yQ(^(iy)  -  ^s(  -  iy))  -  ^^(^(iy)  -  -  iy)) 

-  f  ^-idPziw)  +  tfe(  -  *y)), 

2  ay 

If0  d 

-  /  Kzi(t,x)fi(t)dt  =  i(t/»i(ix)  -  V»i(  -  ix))  +xi— (^i(ix)  -  t/h(  ~  is)),  (49e) 

7r  jd  ax 

-  [  #32 (t,  x)f2(t)dt  =  2 (ifoiix)  +  i)2{  -  ix))  +  x^-(ip2(ix)  +  ip2(  -  ix)).  (49f) 

7r  jd  ax 


(49d) 


The  asymptotic  behavior  of  the  expressions  given  above  near  y  =  0  and  x  =  0  can  now 
be  obtained  by  using  equations  (45)  as  follows: 
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-  f  Ku(t,y)h(t)dt  2  -  2°2_  +  4tt,+i(  _  d)-n Gl(0)(  _  j,)«, 
KJd 


sin(7ro;) 


If" K »(*•  *)*«*  a  -  2Q^a)+1(  -  rf),2G2(°)(  -  »>' 


(50a) 

(50b) 


1  /“u„  /.  \  a  / 1 \  » .  / r)cos('Ka/2)(a  +  2)  +  sin(7ro:/2)(a  +  1) 


7T. 


Ki3(t,y)f2(t)dt  s 


sin(7ra) 


x  (-  a^G^O) (-  y)a, 


(50c) 


l-J°K^,y)h(t)dt  «  (^»(W2)(<»+^-acos(W2)j (  _  a)0G3(0)(  _  y). 


i  f°K3l(t,x)Mt)dt  =s  2 +  i)(  _  d)7,Gi(0)(  -  xf, 

TTJd  sin(7TQ;) 

i  f  K32(t,x)h(t)dt  2  2C0S(J°-^(a  +  2)(  -  4»G2(0)(  -  x)«. 

nJd  Sin(7ra) 


(50d) 

(50e) 

(50f) 


Using  (46)  and  (50),  we  can  express  the  singular  behavior  of  the  integral  equations  (22) 
near  the  end  points  as 

axx(0,y)  =  —  cot(7T7i)(  —  d)aG\{d){y  —  d)71  +  cot(yra)(  -  d)7l6?i(0)(  -  y)a 
2a2  +  4a  4- 1 


sin(7ro:) 


i-dy^m-vT 


+  ( ?7COs(7ra/2)(a  +  2)  +  sin(7ra/2)(a  +  1)  \  ^  _  a^G^0)(  _  y)«; 

\  sin(Tra)  ) 

axy{0,y)  Si  -  cot(7T72)(  -  d)aG2{d)(y  -  d)72  +  cot(7ra)(  -  d)72G2(0)(  -  y)c 
2  a2  +  4a  +  1 


(51a) 


sin(7TO!) 


:(-drG2(0)(-y)a 


+  ( ??sin(7ra/2)(Q  +  1)  -  o:cos(7ra:/2)\  _  _  y ^ 

\  sin(7ro:)  J 


(51b) 
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—  ^  _  a)aG3(a)(x  —  a)^  +  cot(7ra)(  —  0)^73  (0)(  —  x)a 

(k  +  1)  ox 

-  ^^-4g3(x)(  -  x)a(x  -  af 

K  +  1 

+  2S'n{™/f(a  +  1)(  -  d^GjfOX  - *)“ 

sin(7ra) 

+  2C°s(^(a  +  2)(  -  <y  G2(0)(  -  x)“.  (51c) 

sin(7ra) 

The  expressions  given  by  (51),  are  bounded  at  all  points  in  the  intervals  a  <  x  <  0  and 
d  <  y  <  0.  Using  this  condition  characteristic  equations  can  be  obtained  to  determine  the 
unknown  exponents.  If  we  multiply  equation  (51a)  by  (y  —  d)-71  and  let  y  — »  d  and  then 
multiply  by  (  -  y)~a  and  let  y  — »  0,  we  obtain  two  equations  in  the  following  form 


cot(7T7i)  =  0,  (52a) 

7 — y- — -  cos(7ra)  —  (2a:2  4-  4a  +  l)l  (  —  o£)7l(?i(0) 
m(7ra)  L  J 

^cos(7ra/2)(a  +  2)  +  sin(7ra/2)(a  +  1)  (  -  a^G^O)  =  0.  (52b) 


sin(7ra) 

1 


sin(7ra) 


Multiplying  equation  (51b)  by  (y  -  d)-72  and  letting  y  ->  d  and  then  multiplying  by 
(  —  y)~a  and  letting  y  — ►  0  we  obtain  the  following  equations 


COt(7T72)  =  0, 
1 

sin(7ra) 

1 


(53a) 


+  - 


sin(7ra)  L 


cos(7ra)  —  (2a2  +  4a  4-  l)J  (  —  d)72Go( 0) 

?7sin(7ra/2)(a  +  1)  —  acos(7ra/2)l  (  —  a)^Ga(O)  =  0. 


(53b) 


Multiplying  (51c)  by  (x  -  a)~p  and  letting  x  —>  a,  then  multiplying  by  (  -  x)~a  and 
letting  x  — *  0,  following  equations  are  obtained 
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,  ft  —  1 

cot  (717?)  =  -17^-p 


(54a) 


2sin(*V2)  _  „G  +  acos(™/2)  +  2)(  _  d)7IG2(0) 

sin(7ra)  sin(7ra) 


+ 


cot(7ra)  —  r] 


ft  —  li 
ft  +  1. 


-  a)*G3(0)  =  0. 


(54b) 


Using  equations  (52a)  and  (53a)  71  and  72  can  be  obtained  as  —  1/2  and  /3  can  be 
obtained  from  (54a).  As  for  a,  (52b),  (53b)  and  (54b)  constitute  an  eigenvalue  problem 
which  yields  a  characteristic  equation  to  determine  a  and  two  additional  equations  that 
relate  Gi(0),  (?2(0)  and  G3(0).  The  eigenvalue  problem  can  be  expressed  in  the 


following  form  using  (52b),  (53b)  and  (54b), 


011(a)  0  o13(a) 

'J  —  dGi(0) 

'o' 

0  022(a)  o23(a) 

V  -  dG2( 0) 

= 

0 

_a3i(a)  o32(a)  a33(a)_ 

(  -  afG3( 0) . 

0 

The  coefficients  (a)  are  given  in  Appendix  A.  Assuming  Gi(0)  7^  0,  G^O)  7^  0  and 
G3(0)  7^  0  and  solving  the  eigenvalue  problem  given  by  (55)  we  obtain  the  following 
characteristic  equation, 

2a2  +  4a  +  1  -  cos(7ra) 

_ _ _ _ _ L  x 

(ft  +  l)sin2(7ra) 

x  (?7(4a2  4-  10a  +  5  +  (ft  -  l)cos(7ra)  +  ft(2a  +  3))  +  (ft  +  l)sin(7ra))  =  0.  (56) 


The  eigenvalue  problem  (55)  yields  two  more  equations  which  relate  Gi(0),  G3( 0)  and 
G3(0)  as  follows: 


Gi(0)v— d  = 


7Cos(7ra/2)(2  +  a)  +  sin(ira/2)(l  +  a)  ^  /nW  ^p 
2a2+4a  +  l-cos(7ra)  s(  )( 


^2(0)  y/  —  d  — 


7sin(7ra/2)(l  +  a)  —  cos(7ra/2)a 
2a2  +  4a  +  1  —  cos(7ra) 


G3(0)(  -  af 


(57a) 

(57b) 
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For  77  ^  0  the  characteristic  equation  (56)  can  further  be  reduced  to  following  form, 

ri(4a2  +  10a  +  5  +  («  —  l)cos(7ra)  +  «(2a  +  3))  +  («  +  l)sin(7ra)  =  0.  (58) 

In  order  to  verify  the  solution  method  and  the  resulting  characteristic  equation,  the 
analysis  to  determine  the  strength  of  singularity  is  also  carried  out  by  considering  a  90° 
elastic  wedge  and  using  Mellin  transformation.  It  is  shown  that  equation  (58)  can  also  be 
obtained  using  Mellin  transforms.  This  analysis  is  described  in  Appendix  B.  Strength  of 
singularity  a  is  a  function  of  k  and  friction  coefficient  77.  Equation  (58)  can  be  solved  to 
determine  a  for  different  values  of  k  and  rj.  (5  is  the  strength  of  singularity  at  the  other 
end  of  the  contact  area,  and  for  the  case  b  <  0  strength  of  singularities  for  the  contact 
stresses  are  /3  and  u.  For  a  flat  stamp,  the  stresses  are  singular  at  both  ends  of  the  contact 
area  and  for  this  case  equation  (58)  is  solved  for  different  values  of  friction  coefficient 
and  for  v  =  0.25.  The  results  for  a,  j3  and  u  are  shown  in  Figure  6.  If  the  friction 
coefficient  is  greater  than  zero  the  tangential  force  is  directed  as  shown  in  Figure  1  and  a 
is  the  strength  of  singularity  at  the  trailing  end  of  the  contact.  Figure  6  shows  that,  if 
rj  >  0,  a  is  real  and  negative  and  it  is  a  strong  function  of  the  friction  coefficient  77.  For 
the  frictionless  case,  a  =  0  and  no  power  singularity  exists,  and  for  higher  values  of 
friction  coefficient  there  is  a  strong  singularity  at  the  trailing  end  of  the  contact,  u  is  the 
strength  of  singularity  at  the  trailing  end  of  the  contact  for  b  <  0  and  a  is  less  than  u  for 
larger  values  of  friction  coefficient.  77  >  0  is  the  practically  important  case,  because  in  the 
other  case,  i.e.,  77  <  0,  crack  faces  are  forced  to  close  by  the  tangential  force.  In  this  case 
no  singularity  exists  for  6  =  0,  instead  there  is  a  complex  zero,  and  its  imaginary  and  real 
parts  are  shown  in  Figure  6.  The  effect  of  the  Poisson's  ratio  on  the  strength  of  singularity, 
a  is  shown  in  Figure  7.  As  can  be  seen,  the  effect  of  the  Poisson's  ration  is  not  very 
significant  especially  if  77  >  0,  but  the  singularity  becomes  stronger  as  v  decreases.  If 
77  =  0,  no  power  singularity  exists  hence  for  this  case  the  integral  equations  have  to  be 
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checked  in  order  to  determine  the  conditions  not  to  have  a  logarithmic  singularity.  In  this 
case  the  unknown  functions  can  be  expressed  in  the  following  form 


d  <  t  <  0, 

(59a) 

Mt)  =  G2(f)(*  -  dyl/2, 

d  <  t  <  0, 

(59b) 

Mt)  =  G,m  -  af. 

a  <  t  <  0. 

(59c) 

The  asymptotic  behavior  of  the  sectionally  holomorphic  functions  near  point  0  is  found  to 

be, 

■01 0)  =  f — -In  (z), 

7 ry  —  d 

Ip2(z)  s  r — -ln(«). 

Try  “  d 

-  (  (  Lln^)- 

tt(  -  af 

(60) 

The  asymptotic  behavior  of  the  Cauchy  integrals  can  then  be  written  as  follows: 


Coefficient  of  friction,  rj 

Figure  6:  Strengths  of  singularity,  a,  f3  and  u,  u  —  0.25. 
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Coefficient  of  friction,  77 

Figure  7:  Effect  of  the  Poisson's  ratio  on  strength  of  singularity,  a. 


fi(t) 
t  -  y 

m 

t-y 

Mt) 
t-y ' 


dt  £= 


dt  = 


dt  = 


-F.(O) 


7T  y  —  d 

F2(  0) 


Try/  —  d 

Fm 

tt(  -  af 


ln(  -  y), 
ln(  -  y). 


in(-y)- 


(61a) 

(61b) 

(61c) 


Using  the  partial  fractions  (48)  and  (60),  integral  equations  (22)  are  examined  and  it  is 
seen  that  logarithmic  terms  cancel  in  equations  (22a)  and  (22b).  However,  the  analysis 
also  shows  that,  in  equation  (22c)  in  order  to  cancel  the  logarithmic  singularity  following 
condition  has  to  be  satisfied, 

G2{0)\Trd=  —  i<?3(0)(  —  a)*.  (62) 
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Hence,  if  tj  =  0,  there  is  no  power  singularity,  but  there  is  a  relationship  between  G2(0) 
and  G3(0)  as  given  by  (62).* 

4.  Numerical  solution  of  the  singular  integral  equations 

In  this  section,  numerical  solution  methods  will  be  developed  to  solve  the  singular 
integral  equations  given  by  (22)  in  order  to  determine  the  stress  intensity  factors  at  the 
crack  tip  and  contact  stresses.  Three  cases  will  be  considered.  First,  we  will  consider  a 
flat  stamp,  with  b  =  0.  The  strength  of  singularity  a ,  determined  in  section  3.2,  will  be 
used  in  this  solution.  In  the  second  case,  a  flat  stamp  with  b  <  0  will  be  considered. 
Finally,  a  numerical  solution  method  will  be  developed  for  a  circular  stamp  with  smooth 
end  points.  The  numerical  solution  methods  are  based  on  the  expansion  of  the  unknown 
functions  in  Jacobi  polynomials.  Using  the  properties  of  Jacobi  polynomials,  Cauchy 
integrals  can  be  evaluated  in  closed  form.  The  integral  equations  will  then  be  reduced  to  a 
system  of  linear  algebraic  equations  which  can  be  solved  to  determine  the  stress  intensity 
factors  and  contact  stresses. 

4.1  Flat  stamp,  6  =  0 

The  geometry  of  the  problem  is  shown  in  Figure  1.  In  this  section  only  the  problem 
for  b  =  0  will  be  considered.  The  methods  of  solution  for  both  problems  b  =  0  and  b  <  0 
are  similar;  but  in  the  case  of  b  =  0,  the  strength  of  singularity  ft  must  be  considered  and 
there  are  compatibility  relations  between  G\ (0),  C?2 (0)  and  6*3(0)  that  must  be  satisfied 


*  Note  that  for  77  =  0 ,  the  90-degree  wedge  x  <  0,  y  <  0  under  boundary  conditions  v(x,  0)  =  0, 
axy(x,0)  =  0,  crxx(0,  y)  =  0,  axy(0 ,y)  =  0  are  equivalent  to  an  elastic  half  plane  x  <  0  loaded  in  such  a 
way  that  y  =  0  is  a  plane  of  symmetry  and  near  x  =  0,  y  =  0  surface  tractions  axx{0,  y),  crxy{0,  y)  are  zero. 
Hence,  at  and  near  x  =  0,  y  —  0  the  stresses  are  bounded. 
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as  shown  in  section  3.2.  In  order  to  solve  the  problem  numerically,  first  the  integral 


equations  (22)  and  equation  (lg)  are  written  in  normalized  form  as  follows 

-  -  [  +  [  Gu(r,si)gi{r)dr+  f  G13(r,si)g3(r)dr  =  0, 

'Kj-lT-Si  7-1  7-1 

-  -  f  +  [  G22{r,s2)g2(r)dr  +  [  G23(r,s2)g3(r)dr  =  0, 

7T  7-1  T  —  S2  J —l  7-1 

Z"1  AC  1  1  f1  $3(7") 

J  G3l(r,  s3)gi(r)dr  +  J  G32(r,  s3)g2{r)dr  -  ~  ~ 


(63a) 

(63b) 


dr  =  0, 


(63c) 

(63d) 


where,  following 


d  d 
y=2Sl  + 


d 


2’ 

d 


y=2S2+2’ 


a  a 

x  =  2S3  +  r 


and, 


transformations  are  used 
(for  equation  22a), 

(for  equation  22b), 
(for  equation  22c), 


/x  o  ,  { d  d\ 

gi(r)=  -p/i(2r+2  } 

52  (r)  =  \T+^)’ 

53M=  -f/a(f  ’•+§)• 


(64a) 

(64b) 

(64c) 

(65a) 

(65b) 

(65c) 


Note  that,  right-hand  side  of  equation  (63c)  is  zero,  since  for  a  flat  stamp  normal 
displacement  beneath  the  stamp  is  constant.  The  normalized  kernels  Gy  are  given  in 
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Appendix  C.  In  order  to  solve  the  problem,  the  unknown  functions  are  expanded  into 
series  of  Jacobi  polynomials  in  the  following  form 


oo 

91 W  =  (1  -  r)-1/2(l 

71=0 

(66  a) 

OO 

92 ( r)  =  (1  -  r)_1/2(  1  +  ry'£iBnPt1M(r), 

71=0 

(66b) 

OO 

rirMl-rrtHrlW. 

(66c) 

71=0 


Pn  is  the  Jacobi  polynomial,  and  the  exponents  a  and  (5  are  given  by  equations  (56)  and 
(54a)  respectively.  Substituting  (66c)  in  equilibrium  equation  (63d)  and  evaluating  the 
integral  we  obtain  the  first  constant  in  the  expansion  of  gz(r)  in  the  following  form, 


Co  =  -  2/0o, 


2a+p+ir(p  +  i)r(<*  +  i) 

°“  T(a+(3  +  2) 


(67a, b) 


where  T  is  the  Gamma  function.  Substituting  (66)  in  (63)  evaluating  the  Cauchy  principal 
value  integrals  in  closed  form,  truncating  all  the  infinite  series  at  n  =  N,  and  using  (67), 
integral  equations  (63)  can  be  written  in  the  following  form, 

N  ■>  <■>  ->  N 


^A„{ri„(o;)F(n  +  1,  -  n  +  - 

n=0  1  Z 


1  3  1  — si 

2’  2 


a; 


)  +  hnn(si)  |  +  ^  ^C'nfei3n(si) 
J  n=l 


=  —  Cohi3o(si),  —  1  <  si  <  1,  (68a) 

3  1  —  s2  •  '  N 

)  -^  ri22n{S2)  f -r 

n=  1 


N  |  3  .  N 

Y^Bn^r2n(oi)¥(n  +  1,  -  n  +  -  -  a;  +  h22n(s2)j  +  J ~^Cnh2Zn{s2 )  = 

71=0 


-  —  Coh22o(s2), 


-  1  <  s2  <  1, 


(68b) 
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n= 0 


N  N 

^A„/l3in(S3)  +  y>W*3)  + 
n=0 


AT 


l-53 


+  J]C„r3„(a)F(n  +  1,  -  n  -  0  -  a;  1  -  /?;  2 


)  = 


71=1 


l-53- 


-  C'0r3o(a)F(l,  -  0  -  or,  1  -  /?;  -1),  -  1  <  53  <  1, 


(68c) 


where  F()  is  the  Hypergeometric  function,  for  the  definition  and  mathematical  properties 
of  the  Hypergeometric  function  see  for  example  Abramowitz  and  Stegun  [3]  and  for  the 
numerical  computation  of  the  Hypergeometric  function  one  may  refer  to  Luke  [4]. 
Functions  hijn(si),(i,j  =1,2,3)  are  given  in  Appendix  C.  The  formulas  used  in  the 
evaluation  of  the  Cauchy  principal  value  integrals  are  given  in  Appendix  D.  Functions 
rin(z  =  1, 2, 3)  are  given  by, 

2(«-i/2)r(_  1/2)r(n  +  a  +  1) 


(a)  =  r2n(o:)  = 


7rr(n  +  a  +  1/2) 


2(a+flr(/?)r  (w  +  a  +  1) 
Zn  a  irT(n  +  a  +  (3  +  1) 


(69a) 

(69b) 


In  the  solution  of  the  flat  stamp  problem  with  6  =  0,  compatibility  conditions  given  by 
(57a),  (57b)  and  (62)  must  also  be  considered.  Substituting  (66)  in  the  compatibility 
conditions  and  truncating  the  infinite  series,  for  r)  ^  0  we  obtain  the  following  additional 
equations 


J2AnPi  “  1)  = 


a\~a  (\ 


71=0 

N 


-(/m/2)  n 


YJC’‘Pn'")(  -  !). 


n= 0 


,,  -  a/  1  \  -(^+1/2)  N 

£i3„Pi-1/2'«,(-1)  =  R(0)(j)  (i) 

71=0  '  71=0 


(70a) 


(70b) 


where. 
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(71a) 


Pi(a)  = 


77COs(7ra/2)(2  +  a)  +  sin(7ra/2)(l  +  a) 
2a2  +  4a  +  1  —  cos(7ra) 


?7sin(7ra/2)(l  +  a)  —  cos(7ra/2)a 
P2  a  2a2  +  4a  +  1  -  cos(7ra) 


If  rj  =  0,  there  is  only  one  additional  condition 

N 


(71b) 


N  i  ,  _  .  —a  /  1  \  -(/3+1/2)  N 

EB»pi"1/2,“)(  - 1)  =  -  7  G)  (5)  -  «■ 

n=0  '  ' 


(72) 


77=0 


Equations  (68)  can  now  be  numerically  solved  together  with  (70)  or  (72)  to  determine  the 
unknown  constants  An,  Bn  and  Cn.  In  the  numerical  solution  the  equations  are  converted 
to  a  linear  system  of  algebraic  equations  by  using  collocation  points.  If  rj  ^  0,  N 
collocation  points  are  used  for  (68a)  and  (68b)  and  considering  (70)  we  have  two  sets  of 
(TV  +  1)  equations.  Since,  Co  is  known,  (TV)  collocation  points  are  used  for  (68c),  hence 
total  number  of  linear  equations  is  (3TV  +  2)  for  (3TV  +  2)  unknowns.  The  following 
roots  of  the  Chebyshev  polynomials  are  selected  as  collocation  points, 

Sli  =  S2i  =  S3;  =  cos(^ — — — ),  1=1,...,  TV. 


(73) 


Note  that  if  77  =  0,  there  is  one  compatibility  equation  and  in  this  case  (TV  +  1) 
collocation  points  are  used  for  equation  (68a)  and  is  in  the  following  form, 

/7r(2  i  —  1) 


Sli  =  COS 


/tt(2z-  m 
\2(TV  +  1)  / 


i  =  1, TV  +  1. 


(74) 


Stress  intensity  factors  at  the  crack  tip 


Modes  I  and  II  stress  intensity  factors  at  the  crack  tip  are  defined  by 

hi  =  lim y/2(d  -  y)<Txx{0,  y) ,  (75a) 

y— >d 
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(75b) 


ku  =  lim-v/2 (d  -  y)crxy{ 0,  y). 

y-+d 


The  dominant  terms,  for  axx(0,  y)  and  axy{ 0,  y)  near  d  can  be  written  as  follows: 
1  f°h(t)dt 


<7xl{0,V)~  t-y  ■ 
m  ,l[°h(t)dt 

aXy{0,y)  =  -  / 

n  Jd 


t-y 


(76a) 

(76b) 


Considering  the  definition  of  the  sectionally  holomorphic  functions  given  by,  (25  a, b)  the 
asymptotic  behavior  of  (76)  near  d  can  be  written  as, 


v)  =  Gi (d)(  -d)“(d-  y)~1/2, 


(77a) 

(77b) 


Note  that  (77)  is  valid  for,  y  <  d  and  (d  -  y)  -»  0.  Using  the  equations  given  above,  and 
(65a, b),  and  (66a,b)  normalized  stress  intensity  factors  can  be  expressed  as  follows, 

N 


^  V  a  71=0 


(78a) 

(78b) 


After  solving  the  system  of  linear  equations  for  An,  Bn  and  Cn,  stress  intensity  factors 
can  be  obtained  using  equations  (78)  and  contact  stresses  can  be  obtained  using  equation 
(66c). 
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4.2  Flat  stamp,  b  <  0. 


The  geometry  of  the  problem  is  shown  in  Figure  1.  Integral  equations  (22)  can  be 
written  in  normalized  form  as  follows: 


-  -  ['  mi(r)dr  +  f1  Mn(r,sl)ml(r)dr+  [' M13(r,Sl)m3(r)dr  =  0, 
nj-i  r  —  Si  J J- 1 

--[  —  +  [  M22{r,s2)m2(r)dr  +  [  M23(r,  s2)m3(r)dr  =  0, 

ttJ-\  r  —  s2  J- 1  7-i 

/I  1 

Mn(r,  sz)mi(r)dr  +  j  Mz2(r,  sz)m2{r)dr  -  -^-jr]mz{sz) 

_  1  f1  rnzij) 

It  J-i  T  — 

J  mz(r)dr  =  —  2. 


-dr  —  0, 


In  this  case  following  transformations  are  used, 

(for  equation  22a), 

(for  equation  22b), 
(for  equation  22c), 


d  d 
V  —  2Si  +  2’ 

d  d 

y=2S2+r 


a  —  b 

a  +  b 

X  =  2  S3  + 

2  ’ 

and, 

,  ,  b  —  a 

,  /d  d 

mi(r)  =  p 

,  v  6  —  a 

,  (d  d 

m2(r)  =  p 

h{2r+2 

b  —  a 

„  fa-b 

ms(r)  =  p 

H  2  r 

(79a) 

(79b) 


(79c) 

(79d) 


(80a) 

(80b) 

(80c) 


(81a) 

(81b) 

(81c) 


38 


The  normalized  kernels  are  given  in  Appendix  C.The  unknown  functions  are 
expanded  into  series  of  Jacobi  polynomials  as  follows: 

m,(r)  =  (1  -  (82a) 

n=0 

oo 

m2(r)  =  (1  -  r)-1/2 (82b) 

n=0 

oo 

m3(r)  =  (1  -  rf(  1  +  r (82c) 

71= 0 

The  exponents  (5  and  u>  are  given  by  equations  (38)  and  (39)  respectively  and  note  that, 
/5  <  0,  u;  <  0,  /3  +  u=  -1.  Substituting  (82c)  in  equilibrium  equation  (79d)  and 
evaluating  the  integral  we  obtain  the  first  constant  in  the  expansion  of  m3(r)  as  follows: 


Co  =  2sin(7r/5)/7r. 


(83) 


Substituting  (82)  in  (79),  evaluating  the  Cauchy  integrals  in  closed  from  and  truncating 
the  infinite  series  at  N,  we  obtain  the  following  equations, 

N 

-  n+  2’  2’  - 2~^)  +Plln(sl)}  +  TCn^M  — 

n=  1 


=  -  CoPno(si),  -  1  <  si  <  1,  (84a) 

N  1  3  1  _  s  1 

^Bn|f2„F(n+  1,  -  n  +  o’  2’  2 ^  +^P22n(S2)|  +  /  ]CnP2Zn(s2)  = 

n=0 

=  —  CoP23o(s2)>  —  1  <  S2  <  1,  (84b) 

^24„P3in(S3)  +  y~]-BnP32n(53)  +  = 

n=0  n=0  '  n=l 

-  1  <  s3  <  1.  (84c) 
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Functions  Pijn(si),(i,j  =  1,2,3)  are  given  in  Appendix  C.  The  formulas  used  in  the 
evaluation  of  the  Cauchy  principal  value  integrals  are  given  in  Appendix  D  and  Cln  is 
given  by 


r(-i/2)r(n  +  i) 

7t\/2r(n  +  1/2) 


(85) 


By  using  the  collocation  technique  equations  (84),  can  now  be  converted  to  a  linear 
system  of  algebraic  equations.  There  are  a  total  of  (3 N  +  2)  unknowns  in  this  problem. 
Hence  in  the  numerical  solution  (N  + 1)  collocation  points  are  used  for  each  of  the 
equations  (84a)  and  (84b)  and  (N)  collocation  points  are  used  in  equation  (84c).  The 
following  roots  of  Chebyshev  polynomials  are  selected  as  the  collocation  points, 


S2i  =  COS 


/7t(2*  —  1) 
\2(iV  +  1) 


) 


S3,  =  COS 


/7r(2t-  1) 
V  2  N 


)■ 


i  =  1, AT  +  1? 
*  =  1  ,-,N. 


(86a) 

(86b) 


Stress  intensity,  factors  at  the  crack  tip 


Definition  of  the  modes  I  and  II  stress  intensity  factors  are  given  by  equations  (75). 
The  dominant  terms,  for  crIX(0,  y)  and  axy( 0,  y)  near  d  are  given  by  equations  (76).  Using 
these  equations,  asymptotic  behavior  of  the  stresses  near  the  crack  tip  can  be  written  as. 


a„(0,V)sFl(d)(d-y)-'/2, 

(87a) 

Vn(0,y)  =  F2(d)(d-y)~112. 

(87b) 

By  using  the  equations  (81a),  (81b),  (82a),  (82b)  and  (87)  normalized  stress  intensity 
factors  at  the  crack  tip  can  then  be  expressed  as  follows: 
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(88a) 


ku\T^d  -d  p(-i/2to)fl\  (88b) 

P  ib-a)^nn  U> 

After  solving  the  linear  system  of  equations,  stress  intensity  factors  can  be  determined 
using  (88)  and  contact  stresses  can  be  obtained  using  (82c). 

4.3  Circular  stamp 

The  geometry  of  the  problem  is  shown  in  Figure  8.  A  homogenous  half-plane  is 
indented  by  a  rigid  stamp  of  circular  profile.  The  radius  of  the  stamp  is  denoted  as  R,  and 
it  is  assumed  that  contact  area  ( b  —  a)  is  much  smaller  then  the  radius  R.  c  denotes  the 
position  of  the  centerline  of  the  stamp.  In  this  problem  the  contact  area  is  not  constant  and 
it  depends  on  the  force  applied  to  the  stamp.  There  are  four  parameters  in  the  contact 
problem  which  are  P,  a,  b  and  c.  Once  two  of  these  parameters  are  specified,  the  other 
two  have  to  be  determined  from  the  solution  of  the  problem. 
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If  we  assume  that  P  and  one  of  the  end  points  for  instance  a  is  known,  then  b  and  c  can 
be  determined  using  an  iterative  solution  method.  In  order  to  avoid  iterations  and  longer 
computation  time,  in  this  study  a  and  b  will  be  specified  and  corresponding  values  of  P 
and  c  will  be  determined  without  iterations.  Considering  Figure  8,  if  contact  area  is  much 
smaller  than  the  radius  R,  derivative  of  the  normal  displacement  v(x,  0)  can  be  written  in 
the  following  form, 

S-v(x,0)  =  ^-5- a<x<b ,  (89) 

ox  K 


then,  the  integral  equations  given  by  (22)  and  the  equilibrium  equation  given  by  (1.9)  can 
be  written  in  normalized  form  as  follows: 


_1  f1w1(r)dr+  jl  Wu(r,  S\)wi(r)dr  +  ['  W13(r,  Sl)w3(r)dr  =  0, 
7T  J r  —  Si  J_i  J- 1 

--[  W2(r)dL  +  f  W22(r,  s2)w2(r)dr  4-  [  W23(r,  s2)w3(r)dr  =  0, 
nj-i  r-s2  J_i  J- 1 

/I  /*!  ^ 

WZi(r,Sz)wx(r)dr  +  J  W32(r,  s3)w2(r)dr  -  — ^7?W3(s3) 

.if1 2^.- 

nJ_iT  -  s  3  K  +  lv  ii 

IF  L^{r)dr  = 


2  (a  —  b  <x-\-b  2c 
~R  H 


P_ 
/ lR 


(90a) 

(90b) 


(90c) 

(90d) 


In  this  case  following  transformations  are  used, 


d  d 

y  —  2 Sl  +  2’ 

(for  equation  22a), 

(91a) 

d  d 

v=2S2  +  r 

(for  equation  22b), 

(91b) 

a  —  b  d  +  6 

*=  2  ss+  2  - 

(for  equation  22c), 

(91c) 
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and, 


.  ,  1  ,  f  d  d 

Wl  r  =  Jl  1  \2T 2  *’ 


Mr)  =  ±h[\r+d- 


/  \  ^  ( a  ~  b  a  +  b 

w3(r)  =  -h^—r+—^ 


(92a) 

(92b) 

(92c) 


The  normalized  kernels  W{j  are  given  in  Appendix  C.  The  unknown  functions  are 
expanded  into  series  of  Jacobi  polynomials  in  the  following  form, 


oo 

Mr)  =  (1  -  r r1/2f>„.P<-1/2'0)M, 

72=0 

(93a) 

OO 

Mr)  =  (1  -  r), 

71=0 

(93b) 

OO 

«0s(r)  =  (1  -  r/(l  +  rfY^CnP^\r). 

71=0 

(93c) 

The  exponents  /5  and  u>  are  given  by  equations  (38)  and  (39)  respectively  and  in  this  case 
P  >  0,  u  >  0  and  f3  +  u  =  l.  Substituting  (93c)  in  equilibrium  equation  (90d)  and 
evaluating  the  integral,  normalized  force  can  be  expressed  in  the  following  form, 


P_ 

fiR 


a  —  b 
2  R 


CqO  o, 


0o 


2uj(1  -  u)tt 
sin(7ru;) 


(94) 


In  the  numerical  solution,  a/R,  b/R  and  d/R  will  be  specified  and  once  the  problem  is 
solved  the  corresponding  value  of  the  normalized  force  will  be  obtained  by  using  the 
equilibrium  equation,  i.e.,  equation  (94).  In  order  to  determine  c/ R  we  will  consider  the 
consistency  condition.  First  consider  the  last  two  terms  on  the  left-hand  side  of  equation 
(90c),  they  can  be  written  in  the  following  form, 
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.  K-  1  ,  .  If1  U>3  (r) 

P(s3)  =  —7-^77^3(53)  +  r  /  - dr. 


k  +  1 


(95) 


7T  r  —  s3 


If  P(s3)  is  divided  by  the  weight  function  and  integrated  from  -  1  to  1,  following  result 
is  obtained 

rl  F(s3) 


Li - 


■ds3  =  0. 


(96) 


(1  —  s3)/J(  1  +  s3)a 

Other  terms  in  equation  (90c)  must  also  satisfy  this  condition,  and  integrating  the  other 
terms  we  obtain  the  following  expression, 

2  /  a,  —  b %  a b  2c  \  7r 

ynH/32„  =  - 

K 

77—0 


“  ~  ~  -  2  (a-b,  .  a  +  6  2c\  7r 

+  ^  (— C1 "  M  E  )  sin(7ra;)  ’ 

n=0  n=0 


(97) 


where. 


i^3in  =  J  {J  W3i(^,-S3)(l  -  0  1/2  1/2,0)  (r)dr }  (1  -  s3)  ^(1  +  S3)  wds3, 


Wr32n  =  J  {J  W32(r,s3)(l  -  r)  1,2  1/2,0)  (r)dr  j  (1  -  s3)  ^(1  +  s3)  uds3 


(98a) 


(98b) 


Once  the  problem  is  solved  for  a  given  a/R  and  b/R,  c/R  can  be  obtained  using 
equation  (97).  Substituting  (93)  in  (90a,c),  evaluating  the  Cauchy  integrals  in  closed 
from,  eliminating  c/R  from  (90c)  using  the  consistency  condition  given  by  (97)  and 
truncating  the  infinite  series  at  N,  we  obtain  the  following  equations: 

N  ,  o  ,  .  N 


5>n{ftnF( 

71=0 


n  +  1,  — n+ 


1  3  1  —  Sj 


^Pn{0nF(n  +  1,  -  n  +  |  +  922n(s2)}  +  £^923^2)  =  0 

71=0  Ti~~~ 


2  2  2 


1  3  1  s2 


)  +9lln(si)|  +  ]^Cn<7i3n(Si)  —  0, 


n= 0 

1  <  Si  <  1, 

N 


71=0 

-  1  <  s2  <  1, 


(99a) 


(99b) 
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Functions  qijn(si),(i,  j  =  1,2,3)  are  given  in  Appendix  C.  The  formulas  used  in  the 
evaluation  of  the  Cauchy  principal  value  integrals  are  given  in  Appendix  D  and  Vtn  is 
given  by  equation  (85).  Equations  (99)  can  now  be  converted  to  a  linear  system  of 
algebraic  equations  using  collocation  points.  There  are  a  total  of  (3 N  +  3)  unknowns  in 
this  problem,  hence  (N  +  1)  collocation  points  are  used  for  each  of  the  equations.  The 
following  roots  of  the  Chebyshev  polynomials  are  selected  as  the  collocation  points, 

Sli  =  S2j=S3j  =  oos(||^|),  i  =  + 1.  (100) 

For  given  values  of  a/R,  b/R  and  d/R  equations  (99)  can  be  solved  and  normalized 
force  P / fiR  and  c/R  can  be  obtained  using  equations  (94)  and  (97)  respectively. 


Stress  intensity  factors  at  the  crack  tit 


Definition  of  the  modes  I  and  II  stress  intensity  factors  are  given  by  equations  (75). 
Using  the  sectionally  holomorphic  functions,  asymptotic  expressions  can  be  derived  for 
the  stresses  and  the  stress  intensity  factors  can  be  shown  to  be  in  the  following  form, 


(101a) 

(101b) 


Contact  stresses  can  be  obtained  using  equation  (93c). 
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5  Results 


In  this  section,  numerical  results  obtained  for  the  three  different  cases  will  be 
presented  and  the  effect  of  the  friction  coefficient  rj ,  stamp  location  and  crack  length  on 
the  stress  intensity  factors  at  the  crack  tip  and  contact  stresses  will  be  examined.  The 
numerical  solution  methods  are  described  in  Section  4.  Computer  programs  are  developed 
for  the  implementation  of  the  numerical  procedures  and  results  are  presented  in  the 
following  sections  for  flat  and  circular  stamps.  Before,  giving  the  results  for  the  stress 
intensity  factors,  we  examine  surface  stresses  in  the  half  plane  in  the  absence  of  a  crack. 
Figures  9  and  10  show  the  contact  stress  cryy(x,  0)  and  in-plane  stress  axx(x,  0)  for  a  half¬ 
plane  loaded  by  a  circular  stamp.  The  expressions  for  the  stresses  in  the  absence  of  a 
crack  are  obtained  in  closed  form  as  shown  in  Appendix  E.  Figure  10  shows  the  tensile 
peak  for  oxx(x,  0)  at  the  trailing  end  of  the  contact  region.  At  the  leading  end  in-plane 
stress  is  compressive.  The  tensile  peak  at  the  trailing  end  of  the  contact  area  is  the  reason 
of  surface  crack  initiation  in  sliding  contact  problems.  In  the  following  sections  we  will 
present  the  results  for  stress  intensity  factors  and  contact  stresses  for  a  half-plane 
containing  a  crack  at  the  trailing  end  of  the  contact  area. 

5.1  Flat  stamp,  b  —  0 

Figures  11  and  12  show  a  comparison  of  the  results  for  normalized  mode  I  and  II 
stress  intensity  factors  with  those  of  Hasebe  et.  al  [5]  for  the  case  of  a  flat  stamp  with 
6  =  0  and  v  =  0.25.  Mode  I  and  II  stress  intensity  factors  are  given  in  intervals, 
d/a  -  0...1  and  a/d  =  1...0.  Hasebe  et.  al  [5]  solved  this  problem  using  the  complex 
stress  function  technique  and  by  mapping  the  half-plane  with  a  crack  into  a  unit  circle.  As 
also  stated  in  the  paper  by  Hasebe  et.  al  [5],  since  the  half-plane  is  mapped  into  a  unit 
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circle  by  using  a  rational  mapping  function,  the  crack  tips  at  y  =  0  and  y  =  d,  have 
roundness  after  mapping.  However,  as  can  be  seen  in  Figures  11  and  12  the  results  agree 
very  well.  In  the  paper  mentioned,  the  signs  of  the  mode  II  stress  intensity  factors  are 
given  as  negative,  but  for  a  flat  stamp  with  6  =  0,  it  is  clear  that  mode  13  stress  intensity 
factors  have  to  be  positive  at  the  crack  tip.  Figures  1 1  and  12  show  that,  for  low  values  of 
friction  coefficient,  i.e.,  for  a  smaller  tangential  force,  mode  I  stress  intensity  factors  are 
negative.  As  the  tangential  force  increases  mode  1  stress  intensity  factors  increase  and 
they  become  positive.  Hence,  if  the  applied  tangential  force  is  not  sufficiently  large  there 
will  be  contact  between  the  crack  faces  and  crack  closure  will  occur.  In  this  case,  crack 
faces  transfer  compressive  normal  stresses  and  the  assumption  of  zero  crack  surface 
tractions  is  not  valid.  Figure  12  shows  that,  for  a  flat  stamp  with  6  =  0,  mode  II  stress 
intensity  factors  are  positive  for  all  values  of  friction  coefficient  which  is  an  expected 
result  and  if  the  applied  tangential  force  is  smaller,  mode  II  stress  intensity  factor 
increases.  During  numerical  solution,  it  was  observed  that  as  a/ d  or  d/a  approaches  zero, 
it  is  becoming  increasingly  difficult  to  obtain  convergent  results,  but  using  relatively 
larger  number  of  Gaussian  points  and  weights  in  the  numerical  integration  of  the  kernels 
of  integral  equations,  convergent  results  could  be  obtained.  In  Figures  13  and  14  a  flat 
stamp  with  6  =  0  is  considered  and  stress  intensity  factors  are  normalized  with  respect  to 
contact  area  length  a,  in  both  of  the  intervals  d/a  =  0...1  and  a/d  =  0...1.  Hence  for  a 
constant  contact  area  length  a,  the  effect  of  crack  length  d  on  the  stress  intensity  factors  is 
obtained.  It  can  be  seen  that,  as  d  — ►  oo,  both  modes  I  and  mode  II  stress  intensity  factors 
go  to  zero.  It  is  known  that  when  a  brittle  material  is  indented  by  an  indenter  which  at  the 
same  time  translates  laterally  across  the  surface  and  applies  frictional  forces,  tensile 
stresses  intensify  at  the  trailing  edge,  and  this  results  in  the  generation  of  partial  cone 
cracks,  (see  for  example  Lawn  [6]).  The  cone  cracks  extend  into  the  material  and  they 
bend  backwards,  i.e.  opposite  to  the  applied  tangential  force  direction.  An  analysis  of  the 
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results  given  in  Figure  1 4  shows  that  this  is  the  case  for  the  problem  under  consideration. 
Mode  II  stress  intensity  factors  are  positive  for  all  values  of  friction  coefficient  and  crack 
length.  If  we  consider  a  small  element  at  the  crack  tip  as  shown  in  Figure  15,  the  crack 
will  bend  backwards  and  it  will  propagate  in  a  direction  opposite  to  the  applied  frictional 
force. 


Figure  15:  Direction  of  crack  growth 

Figures  16  and  17  show  the  effect  of  the  Poisson's  ratio  on  mode  I  and  II  stress  intensity 
factors  for  a  constant  friction  coefficient  rj  =  0.8.  The  effect  of  the  Poisson's  ratio  on 
mode  I  stress  intensity  factors  is  not  significant  as  shown  in  Figure  16.  The  effect  on 
mode  II  stress  intensity  factors,  however,  is  larger  especially  for  shorter  crack  lengths  and 
as  the  Poisson's  ratio  decreases  mode  II  stress  intensity  factors  increase.  Figures  18  and  19 
show  the  normalized  contact  stresses  for  a  flat  stamp  with  6  =  0.  The  contact  stresses  are 
infinite  at  the  ends  of  the  contact  area  except  for  the  case  77  =  0,  for  which  there  is  no 
singularity  at  the  end  x  =  0. 

5.2  Flat  stamp,  b  <  0 

The  effect  of  the  stamp  location  on  mode  I  and  mode  II  stress  intensity  factors  for  a 
flat  stamp  are  shown  in  Figures  20  and  21.  In  these  figures  (a  —  b)/d  is  kept  constant  as 
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0.1  and  v  =  0.25,  and  stress  intensity  factors  are  given  for  different  values  of  b/ d.  Since, 
stress  intensity  factors  are  normalized  with  respect  to  d,  the  effect  of  the  stamp  location  b 
on  stress  intensity  factors  can  be  clearly  observed  in  these  figures.  The  full  circles  at 
b/d  =  0  are  obtained  using  the  strength  of  singularity  a  and  the  numerical  solution 
method  described  in  section  4.1,  while  the  full  lines  are  obtained  using  the  numerical 
solution  method  described  in  section  4.2.  It  can  be  seen  that,  full  lines  approach  the 
circles  at  b/d  =  0,  which  is  a  verification  of  both  of  the  numerical  methods  developed. 
As  b  increases  from  0,  mode  I  stress  intensity  factors  initially  increase  and  they  go 
through  a  peak  for  larger  values  of  friction  coefficient  rj.  Further  increase  in  b  results  in  a 
decline  in  mode  I  stress  intensity  factors  and  they  approach  zero  for  larger  values  of  b/d 
which  is  an  expected  result.  Similarly,  mode  II  stress  intensity  factors  approach  zero,  for 
larger  values  of  b/d.  It  can  also  be  seen  that,  if  there  is  no  tangential  force,  i.e.,  rj  =  0, 
mode  I  stress  intensity  factors  are  negative  and  mode  II  stress  intensity  factors  are  positive 
for  all  values  of  b/d.  In  Figures  22  and  23  the  effect  of  the  stamp  location  on  mode  I  and 
mode  II  stress  intensity  factors  are  shown  for  (o  —  b)/d  =  1.0  and  u  =  0.25.  Figures  24 
and  25,  show  the  effect  of  the  stamp  location  on  stress  intensity  factors  for 
(a  —  b)/d  =  10.0  and  v  =  0.25.  The  trends  are  similar  and  self-explanatory  in  Figures 
22-25.  Figures  26  and  27  show  the  effect  of  the  Poisson's  ration  on  mode  I  and  mode  II 
stress  intensity  factors  for  (a  —  b)/d  =  1.0  and  rj  =  0.8.  Figures  show  that  for  this  value 
of  (a  —  b)/d,  the  effect  of  the  Poisson's  ratio  is  not  significant,  however  for  a  larger 
value  of  (a  —  b) /d  which  is  taken  as  10.0  in  Figures  28  and  29,  the  effect  of  the  Poisson's 
ratio  is  larger.  Mode  I  stress  intensity  factors  increase  as  v  decreases  as  shown  in  Figure 
28.  The  effect  on  mode  II  stress  intensity  factors  is  shown  in  Figure  29.  Contact  stresses 
for  different  values  of  (a  —  b)/d  are  shown  in  Figures  30-32.  Contact  stresses  are  singular 
at  the  end  points  of  the  contact  area.  Figures  30-32  show  the  effect  of  the  friction 
coefficient  on  contact  stresses 
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5.3  Circular  stamp 


Mode  I  and  II  stress  intensity  factors  for  a  circular  stamp  are  shown  in  Figures  33  and 
34  for  different  values  of  b/R.  In  these  figures  (b  -  a)/R  is  kept  constant  as  0.05  and 
normalized  crack  length  d/R  is  -  0.1.  It  can  be  seen  that  the  variation  of  the  normalized 
stress  intensity  factors  is  similar  to  that  of  the  flat  stamp  case.  Stress  intensity  factors 
approach  zero,  as  b/R  increases.  In  the  solution  of  the  circular  stamp  problem  the  stress 
intensity  factors  are  obtained  for  known  values  of  a/R,  b/R  and  d/R  then  the 
corresponding  value  of  the  normalized  force  P/pR  is  determined  using  the  equilibrium 
equation.  Variation  of  the  normalized  force  with  b/R,  for  (b  —  a)/R  —  0.05  and 
d/R  —  -  0.1  is  shown  in  Figure  35.  It  can  be  seen  that  normalized  force  approaches  a 
constant  value  as  b/R  decreases  and  this  constant  value  can  be  obtained  by  solving  the 
contact  problem  in  a  homogenous  half-plane  without  a  crack,  which  has  a  closed  form 
solution.  This  solution  is  shown  by  the  full  circles  in  Figure  35,  and  our  numerical  results 
approach  the  closed  form  solution  for  small  values  of  b/R.  Closed  form  solution  for  the 
contact  problem  is  given  in  Appendix  E.  As  b/R  goes  to  zero,  normalized  force  increases 
for  smaller  values  of  friction  coefficient  and  it  decreases  for  larger  values.  The  increase 
for  smaller  values  is  due  to  more  compressive  stresses  at  the  crack  tip,  as  the  tangential 
force  decreases  crack  closure  occurs  and  the  normal  force  required  for  contact  increases. 
For  larger  values  of  friction  coefficient,  mode  I  stress  intensity  factors  are  positive  and 
normalized  force  decreases  significantly  with  the  increase  in  b/R.  It  can  also  be  seen  that, 
for  a  given  contact  length  the  normal  force  required  is  smaller  for  a  larger  friction 
coefficient.  Figure  36  shows  the  contact  stresses  for  a/R=  —  0.1,  b/R=  —  0.15  and 
d/R  —  —  0.1.  The  contact  stresses  are  zero  at  the  ends  of  the  contact  area  and  the  effect 
of  the  friction  coefficient  on  contact  stresses  is  shown  in  this  figure.  Another  set  of  results 
for  stress  intensity  factors,  normalized  force  and  contact  stresses  are  shown  in  Figures  37- 
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40,  for  ( b  -  a)/ R  -  0.1.  The  trends  are  similar  to  those  given  in  Figures  33-35.  The 
effect  of  Poisson's  ratio  on  stress  intensity  factors  and  normalized  force  is  also  examined. 
Figure  41  and  42  show  the  modes  I  and  II  stress  intensity  factors  for  different  values  of 
Poisson's  ratio  and  for  q  =  0.8  and  ( b  -a)/R  =  0.1.  It  can  be  seen  that  the  results  are 
very  close  for  all  three  values  of  the  Poisson's  ratio.  However,  Figure  43  shows  that  the 
effect  of  the  Poisson's  ratio  on  normalized  force  is  significant  and  the  normal  force 
required  is  larger  for  a  larger  Poisson's  ratio.  Figures  44-46  show  the  effect  of  the  crack 
length  d/R  on  stress  intensity  factors  and  normalized  force  for  a  circular  stamp  whose 
location  is  given  by  b/R  =  -  0.05  and  a/R  =  -  0.1.  Both  normalized  modes  I  and  II 
stress  intensity  factors  approach  zero  as  d/R  go  to  zero  and  after  a  steep  variation  in 
stress  intensity  factors  for  large  values  of  d/R,  they  eventually  go  to  zero  for  small 
values  of  d/R  as  seen  in  Figures  44  and  45.  Figure  46  shows  the  variation  of  the 
normalized  force  with  crack  length,  it  can  be  seen  that  for  small  friction  coefficients, 
P / \xR  initially  increases  as  crack  length  increases  and  then  goes  to  a  constant  value  as 
crack  length  increases  further.  For  higher  values  of  friction  coefficient  normalized  force 
decreases  initially  and  then  approaches  a  constant  value.  Figures  47-49  show  the 
variation  of  stress  intensity  factors  and  normalized  force  with  crack  length,  for 
b/R  =  -  0.1  and  a/R  =  -  0.05.  The  results  are  similar  to  those  obtained  in  Figures 
44-46. 
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2 xj (o  -  b)  —  (a  +  b)/ (a  —  b ) 

Figure  9:  Contact  stress,  ayy(x,0)  for  a  homogeneous  half-plane  indented  by  a  circular 
stamp  as  shown  in  Figure  8.  v  =  0.25,  d  =  0,  (6  —  a)/f?  =  0.1. 


Figure  10:  In-plane  stress,  &xx(x ,  0)  for  a  homogeneous  half-plane  indented  by  a  circular 
stamp  as  shown  in  Figure  8.  u  =  0.25,  d  =  0,  (6  —  a)/ i?  =  0.1. 
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Figure  11:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  homogeneous  half-plane 
indented  by  a  flat  stamp  as  shown  in  Figure  1,  b  =  0,  v  =  0.25. 


Figure  12:  Mode  II  stress  intensity  factors  for  an  edge  crack  in  a  homogeneous  half-plane 
indented  by  a  flat  stamp  as  shown  in  Figure  1 ,  6  =  0,  v  —  0.25. 
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Figure  16:  Effect  of  the  Poisson's  ratio  on  Mode  I  stress  intensity  factors  for  an  edge 
crack  in  a  homogeneous  half-plane  indented  by  a  flat  stamp  as  shown  in  Figure  1,6  =  0, 
rj  =  0.8. 


Figure  17:  Effect  of  the  Poisson's  ratio  on  Mode  II  stress  intensity  factors  for  an  edge 
crack  in  a  homogeneous  half-plane  indented  by  a  flat  stamp  as  shown  in  Figure  1,6  =  0, 
rj  =  0.8. 
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Figure  18:  Contact  stresses  for  a  homogeneous  half-plane  with  an  edge  crack  indented  by 
a  flat  stamp  as  shown  in  Figure  1,6  =  0,  v  =  0.25,  d/a  =  0.4. 
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Figure  19:  Contact  stresses  for  a  homogeneous  half-plane  with  an  edge  crack  indented  by 
a  flat  stamp  as  shown  in  Figure  1.  v  =  0.25,  a/d  =  0.4. 
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Figure  20:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  homogeneous  half-plane 
indented  by  a  flat  stamp  as  shown  in  Figure  1.  (a  —  b)/d  =  0.1,  v  =  0.25. 


Figure  21:  Mode  II  stress  intensity  factors  for  an  edge  crack  in  a  homogeneous  half-plane 
indented  by  a  flat  stamp  as  shown  in  Figure  1 .  (a  —  b)/d  =  0.1,  v  =  0.25. 


58 


Figure  22:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  homogeneous  half-plane 
indented  by  a  flat  stamp  as  shown  in  Figure  1 .  (a  —  b)/d  =  1 .0,  v  =  0.25. 


Figure  23:  Mode  13  stress  intensity  factors  for  an  edge  crack  in  a  homogeneous  half-plane 
indented  by  a  flat  stamp  as  shown  in  Figure  1 .  (a  — b)/d  =  1.0,  v  =  0.25. 
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Figure  24:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  homogeneous  half-plane 
indented  by  a  flat  stamp  as  shown  in  Figure  1  ,(a  —  b)/d  —  10.0,  u  =  0.25. 


Figure  25:  Mode  II  stress  intensity  factors  for  an  edge  crack  in  a  homogeneous  half-plane 
indented  by  a  flat  stamp  as  shown  in  Figure  1.  (a  -b)/d=  10.0,  u  —  0.25. 
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Figure  26:  Effect  of  the  Poisson's  ratio  on  mode  I  stress  intensity  factors  for  an  edge 
crack  in  a  homogeneous  half-plane  indented  by  a  flat  stamp  as  shown  in  Figure  1. 
(a  —  6)/d  =  1.0,  77  =  0.8. 


Figure  27:  Effect  of  the  Poisson's  ratio  on  mode  II  stress  intensity  factors  for  an  edge 
crack  in  a  homogeneous  half-plane  indented  by  a  flat  stamp  as  shown  in  Figure  1. 
(a  —  b)/d  =  1.0, 77  =  0.8. 
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Figure  28:  Effect  of  the  Poisson's  ratio  on  mode  I  stress  intensity  factors  for  an  edge 
crack  in  a  homogeneous  half-plane  indented  by  a  flat  stamp  as  shown  in  Figure  1. 
(a  — b)/d  =  10.0 ,77  =  0.8. 


Figure  29:  Effect  of  the  Poisson's  ratio  on  mode  H  stress  intensity  factors  for  an  edge 
crack  in  a  homogeneous  half-plane  indented  by  a  flat  stamp  as  shown  in  Figure  1. 
(a  —  b)/d  =  10.0,  rj  =  0.8. 
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a  flat  stamp  as  shown  in  Figure  1.  v  =  0.25,  b/d  =  0.4,  (a  —  b)/d  =  1.0 
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Figure  32:  Contact  stresses  for  a  homogeneous  half-plane  with  an  edge  crack  indented  by 
a  flat  stamp  as  shown  in  Figure  1.  v  —  0.25,  b/d  —  0.4,  (a  —  b)/d  =  10.0. 
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Figure  33:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  homogeneous  half-plane 
indented  by  a  circular  stamp  as  shown  in  Figure  8.  (b  —  a)/R  =  0.05,  d/R  =  —  0.1, 
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Figure  34:  Mode  II  stress  intensity  factors  for  an  edge  crack  in  a  homogeneous  half-plane 
indented  by  a  circular  stamp  as  shown  in  Figure  8.  (b  -  a)/ R  =  0.05,  dJR—  -  0.1, 
v  =  0.25. 


b/R 

Figure  35:  Normalized  force  for  a  homogeneous  half-plane  containing  an  edge  crack  and 
indented  by  a  circular  stamp  as  shown  in  Figure  8.  (b  —  a)/R  =  0.05,  d/R  —  —  0.1, 
v  =  0.25. 
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Figure  36:  Contact  stresses  for  a  homogeneous  half-plane  with  an  edge  crack  indented  by 
a  circular  stamp  as  shown  in  Figure  8.  v  =  0.25,  d/R=  —0.1,  b/R=  —0.1, 
a/R  =  -0.15 


Figure  37:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  homogeneous  half-plane 
indented  by  a  circular  stamp  as  shown  in  Figure  8.  (b—  a)/R  =  0.1,  d/R  =  —  0.1, 
v  =  0.25. 
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b/R 

Figure  38:  Mode  II  stress  intensity  factors  for  an  edge  crack  in  a  homogeneous  half-plane 
indented  by  a  circular  stamp  as  shown  in  Figure  8.  (6  —  a)/R  =  0.1,  d/R=  -  0.1, 
v  =  0.25. 


b/R 

Figure  39:  Normalized  force  for  a  homogeneous  half-plane  containing  an  edge  crack  and 
indented  by  a  circular  stamp  as  shown  in  Figure  8.  (b  -  a)/ R  =  0.1,  d/R  =  -  0.1, 
v  =  0.25. 
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Figure  40:  Contact  stresses  for  a  homogeneous  half-plane  with  an  edge  crack  indented  by 
a  circular  stamp  as  shown  in  Figure  8.  v  —  0.25,  d/R=  —0.1,  b/R=  —0.1, 
a/R  =  -  0.2. 


Figure  41:  Effect  of  the  Poisson's  ratio  on  mode  I  stress  intensity  factors  for  an  edge 
crack  in  a  homogeneous  half-plane  indented  by  a  circular  stamp  as  shown  in  Figure  8. 
(b  -a)/R  =  0.1,  rj  =  0.8,  d/R  =  -0.1. 
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Figure  42:  Effect  of  the  Poisson's  ratio  on  mode  II  stress  intensity  factors  for  an  edge 
crack  in  a  homogeneous  half-plane  indented  by  a  circular  stamp  as  shown  in  Figure  8. 
(b  -a)/R  =  0.1,  rj  =  0.8,  d/R  =  -  0.1. 


b/R 

Figure  43:  Effect  of  the  Poisson's  ration  on  the  normalized  force  for  a  homogeneous  half¬ 
plane  containing  an  edge  crack  and  indented  by  a  circular  stamp  as  shown  in  Figure  8. 
(i b  -  a)/R  =  0.1,  d/R  =  -  0.1,  rj  =  0.8. 
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Figure  44:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  homogeneous  half-plane 
indented  by  a  circular  stamp  as  shown  in  Figure  8.  b/R  —  -0.05,  a/R—  -0.1, 
v  =  0.25.  " 


Figure  45:  Mode  II  stress  intensity  factors  for  an  edge  crack  in  a  homogeneous  half-plane 
indented  by  a  circular  stamp  as  shown  in  Figure  8.  b/R  —  —0.05,  a/R=  —0.1, 
v  —  0.25.  " 
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Figure  46:  Normalized  force  for  a  homogeneous  half-plane  containing  an  edge  crack  and 
indented  by  a  circular  stamp  as  shown  in  Figure  8.  b/R=  —0.05,  a/R=  —0.1, 
v  =  0.25.  * 


Figure  47:  Mode  I  stress  intensity  factors  for  an  edge  crack  in  a  homogeneous  half-plane 
indented  by  a  circular  stamp  as  shown  in  Figure  8.  b/R=  —0.1,  a/R=  —0.15, 
v  =  0.25.  " 
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Figure  48:  Mode  II  stress  intensity  factors  for  an  edge  crack  in  a  homogeneous  half-plane 
indented  by  a  circular  stamp  as  shown  in  Figure  8.  b/R=  —0.1,  afR—  —0.15, 
v  —  0.25.  ’ 
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Figure  49:  Normalized  force  for  a  homogeneous  half-plane  containing  an  edge  crack  and 
indented  by  a  circular  stamp  as  shown  in  Figure  8.  b/R=  —0.1,  a/R—  —0.15, 
v  =  0.25.  ' 
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APPENDIX  A 


Coefficients  for  the  eigenvalue  problem 


Following  are  the  coefficients  for  the  eigenvalue  problem  to  determine  the  strength  of 
singularity  a  (see  equation  (55)). 


an  (a)  =  a22(a)  =  — - — -fcos(7ra)  -  (2a2  +  4a  +  l)), 
sin(7ro;)  V  / 

(Ala) 

ai3(a)  =  sin^7rQj  ^cos(7ra/2)(2  +  a)  +  sin(7ra)(l  +  a)) , 

(Alb) 

«23(«)  =  (^sin(7ra/2)(l  +  a)  -  acos(7ra/2)), 

(Ale) 

2sin(7ra/2)(l  +  a) 

^31  (^J  —  •  /  \  ? 

sm(7ra) 

(Aid) 

.  2cos(7ra/2)(2  +  a) 

&32  (^)  —  •  /  \  ? 

sm(7ra) 

(Ale) 

.  .  ,  .  K  -  1 

a33(a)  =  cot(Tra)  rj. 

K  “T"  1 

(A  If) 

73 


APPENDIX  B 


Mellin  transform  method  for  the  derivation  of  the  characteristic  equation 

The  characteristic  equation  given  by  (58)  can  also  be  obtained  by  considering  a  90- 
degree  elastic  wedge  as  shown  in  Figure  50  and  using  Mellin  transformation. 


Figure  50:  90-degree  elastic  wedge 
The  boundary  conditions  of  the  problem  are  given  by, 


are(r,0)  -rjaee{r,0)  =  0, 

0  <  r  <  oo, 

(Bla) 

-j^ue(r,0)  =  /(r), 

0  <  r  <  oo, 

(Bib) 

ag6(r,n/2)  =  0, 

0  <  r  <  oo, 

(Blc) 

■? 

Qi 

to 

II 

o 

0  <  r  <  oo. 

(Bid) 

In  polar  coordinates,  following  definition  of  the  stresses 

identically  satisfies  the  equilibrium  equations, 

in  terms  of  a  stress  function  x 

a  Ir  e)-19x(r’e)  +  152*(r’ 

arr{T,V)  r  dr  1  r2  QQ2 

0) 

5 

(B2a) 

r_  n)  ...  d2x(r,  9) 
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(B2b) 

Thus,  the  compatibility  condition  becomes. 


V2V2xM)  =0, 


(B3) 


where, 


d2  id  id2 


v  =  — | - h 


(B4) 


dr2  r  dr  r 2  d92 

It  is  also  known  that  (see  for  example  Hein  and  Erdogan  [7]),  displacements  can  be 
expressed  in  terms  of  this  biharmonic  stress  function  x  and  another  harmonic  function  4> 
in  the  following  form, 


dx(r,6)  ,  /«  + 1\  d<p{r,6) 


o  /  /v\  ldx(r,0)  ,  /«+l\  2 d<t>{r,6) 

2W(r,»)=  ---5r-  +  (— jr-^, 

x(r,  9)  and  </>(r,  6)  are  related  by 

and  0(r,  6)  is  a  harmonic  function 
V24>(r,  9)  =  0. 

Mellin  transform  of  a  function  /(r)  and  its  inverse  are  defined  as 

^  1  pc-\-ioo^ 

7(P )  =  /  f{r)rv~ldr,  /(r)  =  —  /  f(p)r~pdp. 

JO  *mJc-ioo 

Mellin  transform  of  the  derivative  can  be  expressed  as  follows: 


(B5a) 

(B5b) 


(B6) 


(B7) 


(B8) 


(B9) 


provided  that 


JTO-i 

rp+m~1^f(r)  =  0, 

as,  r  — *  0  and  r  — »  oo  for  m  =  1,  2, 3, n. 


(BIO) 


Taking  the  Mellin  transform  of  equations  (B3),  (B6)  and  (B7),  following  equations  are 
obtained, 


dAx(0,P)  ,  (_*  ,  ,  on2 \d2x(Q,p)  ,  „2/_  ,  0X2 


+  {p2  +  (p  +  2)2)- 


dQA  '  '  v-  '  "/  /  d62 

+  p2<t>(d,p)  =  0, 


+  p  (p  +  2)  x(0,p)  =  0, 


d24>{0,p)  ,  2 


de 2 


(Blla) 

(Bllb) 

(Bile) 


Solving  equations  (Blla)  and  (Bllb)  and  using  (Bile)  displacements  and  stresses  can  be 
expressed  in  the  following  form, 


X 


dur  .  due ' 


Oue\  _ 
dr  '  dr  J 

i  pc+io o 

-  —  /  (p  +  1)  (Apeip6  +  B(p  +  l)ei(p+2)(?  +  KBe~i{j,+2)e)r-pdp, 
2  ™Jc-ioo 


(B12a) 


l 

^7  (o’rO  +  ZCrfltf)  = 


2 h£+lC°(p  +  ^  ( ApeiP$  +  +  iyiP+2)e  -  Be~i(p+2)ey-pdp, 


(B12b) 


In  equations  (B12),  A  and  B  are  complex  constants  and  B  is  the  complex  conjugate  of  B. 
Using  equations  (B12)  and  boundary  conditions  (Bl)  following  expressions  are  obtained 
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o«(r,°)  _  r  1  r+i“  ^  +  to  +  1  -  °~(™)  S_da< 


4^  jo 
i  d 


2m. 


D(a) 


ta+l 


(B13a) 


2(«  +  1)  <9r 


u0(r,ir/2) 


poo  1  pc-t 

l  mdt^L 


1  fc+lco  rjcos(na)(a  +  2)  +  sin(7ra/2)(a  +  1)  ra 


L»(a) 


ta+1 


da,  (B13b) 


where, 

D(a)  =  77(4a2  +  10a  +  5  +  (k  -  l)cos(7ra)  4-  «(2a  +  3))  +  («  +  l)sin(7ra).  (B14) 

Note  that  £>(a)  is  same  as  the  characteristic  equation  (58).  If  one  performs  the  inversions 
in  (B13)  by  using  the  theory  of  residues,  then  negative  roots  of  D(  a)  give  stresses  that  are 
singular  as  r  approaches  zero.  Thus,  D{a)  =  0  is  the  characteristic  equation  of  the 
problem. 
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APPENDIX  C 


Kernels  of  the  integral  equations 

The  integral  equations  for  are  given  by  and  the  kernels  are  given  by  (22)  and  (23) 
respectively.  In  this  section  we  will  give  the  transformed  form  of  the  kernels  and  other 
terms  for  the  three  cases  considered. 


Flat  stamp,  b  =  0 

The  transformed  form  of  the  kernels  used  in  equations  (63  a, c)  are  given  by, 


£n(r,  Si)  = 

1  cf  , 

?d  d  d  d' 

—  7T  2Kn 

.2r+2’2S,  +  2. 

Gu(r,si)  = 

1  a  T,  / 

f  a  ad  d 

-^“1 

,2r+2'2Sl  +  2 

G22{r,s2 )  = 

1  d  T,  1 

( d  d  d  d 

^2r+2’2S2+2 

£23(^5  s2)  = 

1  a  „  i 

'a  ad  dy 

-*2K^ 

.2r+2’2S2+2- 

£31  (^5  S3)  = 

1  d„  , 

^d  d  a  a~ 

,5r+2’2S3+2. 

£32(^5  s3)  = 

Id  „  / 

f d  da  a' 

,2r+2'2S3+2. 

The  terms  used  in  equations  (68)  are  in  the  following  form 
hnnisx)  =  £cn(r,Sl)(l  -  r^l  +  r)aPt1/2’a\r)dr, 

hl3n(si)  =  f_G13(r,Sl)(  1  -  rf(l  +  r)aP^(r)dr, 

^22n(s2)  =  J  £22(7-5  S2XI  ~  r)  ^2(1  +  r)ap(  l!2,a\r)dr, 


(Cla) 

(Clb) 

(Clc) 

(Cld) 

(Cle) 

(Clf) 

(C2a) 

(C2b) 

(C2c) 
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ha n(s2)  =  J'  cai r, »,)( 1  -  +0  +  r)‘P*f^(r)dr,  (C2d) 

/*,„(«»)  =  f  G3l(r,s3)(l  -  rrl'2(l  +  rrPt'M(r)dr,  (C2e) 

*>*,(*>)  =  fG^sMl-ry^il  +  rTPl-U^dr.  (C2t) 

Flat  stamp,  b  <  0 

The  transformed  form  of  the  kernels  used  in  equations  (79a, c)  are  given  by, 

.  .  1  d  T_  ( d  d  d  d\  { r^,~i r,\ 

Mn(r,s,)=  ---Kll{-r+-,-s,  +  ~),  <C3a) 

1  a  —  b  (a  —  b  a  +  b  d  d\  /oou\ 

M13(r,Sl)=  ~-—K^y—r+—,-s1  +  -J,  (C3b) 

.  .  I  d  ( d  d  d  d\  x 

M22(r,s2)=  -;2if22(2r+2’2S2+2j’  (C3C) 

,  1  a-b  /a-b  a  +  6  d  d\ 

M23(r,52)=  ^2s(— r+— ,-52  +  -],  (C3d) 

, ,  .  ,  1  dT.  (d  d  a-b  a  +  b\  fn-  . 

Af3i(r,s3)=  --2Jf!1(2r+2’_2-S3  °3 

. ,  .  .  1  dr.  (d  d  a-b  a  +  b\ 

M32(r,s3)=  --2^32(2r+2’_T-S3  +  ^_)-  (C  ^ 

The  terms  used  in  equations  (84)  are  in  the  following  form 

fti»(si)  =  J  Mn(r,  si)(l  -  ry1/2Pt1/2fi)(r)dr ,  (C4a) 

Pi3n(si)  =  f_Ml3(r,  5l)(l  -  r)*(l  +  rfP^\r)dr,  (C4b) 

P22n(S2)  =  JlM22(r ,  S2)(l  -  r)-1/2p(-1/2.°)(r)dr.  (C4c) 
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JtotaM  =  «2)(1  -  r)"(l  +  r  fP^\r)dr,  (C4d) 

J*l»(<3)  =  £m3i( r,<.8)(l  -  r)-l/2Pt1/2'0)(r)dr,  (C4e) 

P32„M  =  £  M32(r,  s3)(l  -  r)-*/»p,(- W)(r)dr.  (C4f) 

C.3  Circular  stamp 

The  transformed  form  of  the  kernels  used  in  equations  (90a, c)  are  given  by, 

rrr  .  .  Id  r.  /d  d  d  d\  N 

W,n(r,si)=  -;2ifll(2r+2’2Sl  +  2)’  <C5a) 

W'i3(r,s,)=  ■ K*[— r  +  —  '2Sl  +  2j’  (C5b> 

S2)  =  -  + -, -S2  + 

„TT  .  .  1  a  —  bT./a  —  b  a  —  bd  d\  ,nfJ, 

W23(r,s2)=  -  -  -  K2z (  2  r  ^ - 2~ ’2S2+2/’  (C5d) 

Tir  .  .  1  d  „  (d  d  a -b  a-b\  \ 

^(r.sa)-  "  ^  5^1  (jr  + — s3  +  — ),  (C5e) 

rxr  /  x  1  d  ts  fd  d  a -b  a-b\ 

W*>(.r,*)  =  —  7r2^32\2r"*”  2’  2  S3  2  )*  (C5° 

The  terms  used  in  equations  (99)  are  in  the  following  form 

«n„(*l)  =  J'w „(r,  Sl)(l  -  r)-1/2pyi22.°>(r)dr,  (C6a) 

«ita(»l)  =  £w18(r,»1)(l-ry,(l  +  r)"i^)(r)<ir>  (C6b) 

qnjs2)  =  j  Wn{r,si){\  -  r)-1/2Pt1/m(r)ir,  (C6c) 
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?23n  (^2)  = 

fw23(r,  s2)(l  -  rf(  1  +  rfP^\r)dr, 

(C6d) 

<?31n(53)  = 

1  W3i(r,s3)(l  -  ryl/2Pjr1/2fi)(r)dr, 

(C6e) 

932ra(53)  = 

1  W32(r,  s3)(l  -  r)_1/2P7|_1/2,0) (r )dr. 

(C6f) 
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APPENDIX  D 


Closed  form  expressions  for  Cauchy  principal  value  integrals 


Following  result  which  is  given  by  Tricomi  [8],  is  the  basic  expression  that  is  used  to 
evaluate  the  Cauchy  principal  value  integrals  involving  Jacobi  Polynomials, 


7T  J  _  i  LX 

cot(tro;)(l  —  rr)a(l  +  xYP^\x)  — 
2(a+/5)r(a)r(n  +  (3+l) 


f(h  +  1.  -  n  -  a  -  /3]  1  -  a;  ^ 


7rr(n  +  a  +  (3  +  1) 


(Dl) 


where  a  >  -  1 ,  (3  >  -  1,  a  ^  0, 1, 2, T  is  the  gamma  function  and  F()  is 
hypergeometric  function.  If  (a  +  (3)  is  equal  to  —1,0  or  1,  (Dl)  can  be  further 
simplified  as  follows. 


■nj- 1  t-x 


2~x 


=  COt(TO)(l  -  *)“(!  +  xfP^(x)  ~  ^X-^pjr^x), 


sin(7ra) 


(D2) 


where,  x  =  —  (a  +  (3). 
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APPENDIX  E 


Indentation  of  a  homogeneous  half-plane  by  a  circular  stamp 

The  expressions  for  the  contact  stress  oyy(x,  0),  in-plane  stress  axx(x,0)  and 
normalized  force  P / fiR  can  be  obtained  in  closed  form  for  the  problem  shown  in  Figure 
51.  The  governing  equations  can  be  written  as, 


Figure  51 :  Homogeneous  half-plane  indented  by  a  circular  stamp 


f 


4  fi  x  —  c 

K  -f-  1  R 


a  <  x  <  b, 


cryy(x,  0 )dx  —  —  P. 


(Ela) 

(Elb) 


(Ela)  and  (Elb)  can  be  solved  in  closed  form  and  following  results  can  be  obtained 


x*  -  2x/{a  —  b)  —  (a  +  b)/ (a  —  b ), 

2sin(7r/3)  (a  -  b)  /  *\0fi  ,  *Y 

°yy(x  »°)  -P  K  +  1  R  (2  x)  (1  +  *) 


-  1  <  x*  <  1, 


(E2a) 

(E2b) 
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+  -  2n  rzb)20(l  0), 

fiR  k  +  1  v  R  ) 

(E2c) 

c  =  |{< 

a-b)(2P-  l)  +  (o  +  6)}, 

(E2c) 

where, 

cot  (7 r/3) 

II 

1 

-3 

*  'ST 

+  1 

1-1 

(E3a) 

COt(7Rt>) 

*?(«  -  !) 

K  +  1 

(E3b) 

and,  P  >  0,  u  >  0,  P  +  u  =  1.  It  can  be  seen  that  if  a  and  b  are  known,  normalized  force 
can  be  obtained  using  (E2b)  and  position  of  the  centerline  of  the  stamp  can  be  obtained 
from  (E2c).  The  in-plane  stress  axx{x,  0)  can  also  be  determined  in  closed  form  (for  the 
details  of  the  derivation  of  uxx(x,  0),  see  Guler  [9]).  For  —  oc  <  x*  <  —  1,  we  can 
write, 

axx(x*,  0)=  -C1(-x*  +  lf( -x*-l)w-x*+p-u,  (E4a) 


for  -  1  <  x*  <  1, 

axx(x *,  0)  =  Ci(l  -  x*f{l  +  x*)wcos(nP)  -  x*  +  p  -  u  +  C0(  1  -  x*f(l  +  x*f , 

(E4b) 


and  for  1  <  x*  <  oo, 


axx(x\ 0)  =  CPx*  -  1)  V  +  1  )u-x'  +  P-u, 


(E4c) 


where 


„  (a  -  6)  2sin(7 r/3) 

Go  =  H  5 

R  K  +  1 


Ci  = 


2Cp 

sin(7r/5) 


(E5a,b) 
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